NPS  ARCHIVE 
1968 
BUSH,  T. 


LARGE  ANGLE   SCATTERING   OF 
LITHIUM  IONS   BY   HELIUM  ATOMS 


J3R&RY 
E4VAL  FQS 

,  GALIF. 


UNITED  STATES 
NAVAL  POSTGRADUATE  SCHOOL 


THESIS 


LARGE  ANGLE   SCATTERING  OF 
LITHIUM   IONS    BY  HELIUM  ATOMS 
by 

Thomas   Oliver    Bush 


September    1968 


This  document  has  bp<=>n  apprcrtd  for  publio 
release  and  sale;   i       ti     trlfcctlon  !s  unlioited- 


LARGE  ANGLE  SCATTERING  OF 
LITHIUM  IONS  BY  HELIUM  ATOMS 

by 


Thomas  Oliver  Bush 

Lieutenant,  United  Spates  Navy 

B.S.,  University  of  Oklahoma,  i960 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

DOCTOR  OF  PHILOSOPHY 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
September  1968 


abstract!  J 
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An  experiment  using  the  focusing  properties  of  a  nonuniform 
axially  symmetric  magnetic  field  to  study  large  angle  (9  >  32°) 
scattering  of  ions  from  atoms  and  molecules  is  described.   The 
scattering  cell  is  placed  on  the  magnetic  field  axis  and  all  ions 
in  the  momentum  interval  I Apl  at  |pj  scattered  into  a  conical  shell 
AO  at  0  are  counted  by  a  detector  placed  further  down  the  axis.  The 
solid  angle  is  increased  by  a  factor  of  several  hundred  over  con- 
ventional scattering  techniques  of  comparable  angular  resolution. 
The  magnetic  vector  potential  and  the  trajectory  equation  for  a 
charged  particle  in  the  field  are  derived.   The  computer  program 
which  solves  the  trajectory  equation  is  given  along  with  a  detailed 
description  of  the  apparatus.   Our  measurements  of  the  absolute 
differential  scattering  cross  sections  for  Li   on  He  from  36   to 
41.7°  at  200  eV  and  300  eV  and  for  He+  on  He  from  40°  to  54°  at 

300  eV  and  400  eV  demonstrate  the  feasibility  of  the  technique. 

+ 
Usxng  the  measured  cross  sections  for  Li   on  He  and  the  computa- 
tional method  of  Firsov  we  have  calculated  the  interatomic  poten- 
tial function  for  nuclear  separations  of  O.15X  to  1.1A*, 
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CHAPTER    I 
INTRODUCTION 
In   recent   years    the   techniques    for    measuring   differential 
scattering    cross    sections    for    various    atomic   and  molecular    pro- 
cesses   have   been    improved  to   the  point   where   very  precise  mea- 
surements   can   be   made    in   an   almost   routine   manner .       In    the   con- 
ventional  approach   a    target    of   N     particles    is    confined   within   a 
surrounding    vacuum.      A  well    collimated  beam   of   nearly   monoener - 
getic    ions    of   current    I.    is   projected   into   the   target.      A   de- 
tector,   subtending   a    solid  angle   dQ  (9)    =    Sin   0   dO   d0   at    the 
interaction   region,    is   placed   at    an   angle   0   to   the    incident    ion 
beam  and   the   number    of  particles    entering    the   detector,.  N    ,    is 
measured.      The   scattered  particles    are  usually   detected   as   a 
current   where   I      =   N   dQ(O).      From   these  measured   quantities,    the 
differential    scattering   cross    section    in    the    laboratory    coordinate 
system   can   be   calculated  as 

da(0)  =1/1  .NdQ(O)  . 
v     '  sit 

However,  in  order  to  obtain  reasonably  good  angular  resolution, 

o  s    o 

the  solid  angle  includes  only  about  1   of  the  available  360   of 

azimuthal  angle  0  into  which  ions  scatter.   This  restriction  re- 

-k 
suits  in  fairly  small  solid  angles  of  about  10    steradians, 

-20  2 
which  sets  a  lower  limit  of  around  10  cm  on  observable  labora- 
tory differential  scattering  cross  sections.  By  modifying  the 
techniques  of  thin  magnetic  lens  p-ray  spectroscopy,  we  have  con- 
structed an  experimental  apparatus  which  collects  essentially  all 
ions  scattered  into  a  360°  cone  dO  wide  at  0.  The  solid  angle  in 
this  apparatus  is  about  300  times  larger  than  that  of  other  methods, 

11 


allowing  us  to  observe  differential  scattering  at  laboratory 
angles  greater  than  32  degrees  and  cross  sections  substantially 
smaller  than  any  previously  measured, 

The  apparatuss  shown  schematically  in  Figure  1,  was  devel- 
oped to  make  use  of  the  focusing  properties  of  an  axially  sym- 
metric nonuniform  magnetic  field.   The  scattering  cell  exit  slit 
is  open  290   of  the  aximuthal  angle,  allowing  ions  of  a  given  p 
and  9  but  scattered  at  most  any  azimuthal  angle  to  be  focused 
at  the  detector  by  the  magnetic  field,   Both  the  magnet  current 
and  the  distance  between  the  scattering  cell  and  detector  are 
variable.   Various  combinations  of  the  magnet  current  and  the 
detector  distance  make  a  wide  range  of  reactions  observable  with 
the  apparatus „   The  details  of  the  apparatus  are  discussed  in 
Chapter  III. 

The  differential  equation  for  the  trajectory  of  a  charged 
particle  in  a  nonuniform  axially  symmetric  magnetic  field  and 
the  magnetic  vector  potential  for  this  field  are  developed  in 
Chapter  II.   The  trajectory  equation  is   numerically  integrated 
to  give  the  values  of  the  magnetic  field  and  detector  position 
necessary  to  observe  charged  particles  scattered  at  various  angles 
and  energies.   The  solid  angle  is  calculated  from  the  computed 
trajectories  and  is  a  function  of  the  angle  of  scatter,  the  en- 
ergy of  the  scattered  ion,  the  magnetic  field,  and  the  detector 

distance.   This  calculation  is  discussed  in  detail  in  Chapter  II. 

+ 
The  elastic  differential  scattering  cross  sections  for  He 

on  He  at  300  eV  from  40°  to  52°,  and  for  He+  on  He  at  400  eV  from 

42   to  53   were  measured  using  this  method  and  the  results  are  in 
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essential  agreement  with  measurements  of  Lorents  and  Aberth  who 
used  a  standard  technique.   The  elastic  differential  scattering 
cross  sections  measured  for   Li   on  He  at  200  eV  and  300  eV  from 

34   to  4l.7   also  compared  favorably  to  extrapolated  values  of 

2 
cross  sections  measured  at  smaller  angles  by  other  investigators  . 

These  experimental  results  are  given  in  Chapter  IV, 

Finally,  the  elastic  differential  scattering  cross  section 

6  + 
for   Li   on  He  measured  by  our  method,  along  with  data  taken  at 

smaller  angles  by  other  workers,  was  used  to  calculate  the  ion- 
atom  potential  function  for  this  system  using  a  method  developed 

3  . 

by  Firsov  .   This  calculation  is  discussed  in  Chapter  V. 

Based  on  these  results  we  conclude  that  the  feasibility  of 

the  technique  has  been  established  and  its  usefulness  in  obtaining 

new  data  has  been  demonstrated. 
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CHAPTER  II 
KINEMATICS  OF  CHARGED  PARTICLES  IN  A  NONUNIFORM 
AXIALLY  SYMMETRIC  MAGNETIC  FIELD 
In  this  chapter  the  equation  of  motion  of  a  charged  particle 
in  a  uniform  magnetic  field  is  solved  to  show  the  explicit  re- 
lation among  the  following  three  variables:  the  vector  momentum 

— »  — » 

of  the  ion,  p,  the  magnetic  field,  B,  and  Z  ,  the  distance  from 

the  origin  at  which  a  detector  must  be  placed  to  observe  the  ion. 

The  expression  for  the  solid   angle  is  obtained  in  closed  form 

and  its  dependence  on  the  various  parameters  is  shown. 

The  orbit  equation  for  the  axially  symmetric  nonuniform 
magnetic  field  is  then  derived  along  with  the  expression  for  the 
magnetic  vector  potential,  which  is  part  of  the  orbit  equation. 
Expressions  for  the  angular  and  energy  resolution  are  found  and 
these  results  are  used  to  find  the  solid  angle.   The  dependence 
of  the  solid  angle  on  the  various  parameters  is  shown  and  dis- 
cussed. 

Finally,  an  experimental  method  of  simulating  a  particle 
trajectory  in  a  magnetic  field  using  a  flexible,  current  carrying 
wire  is  described.   The  results  of  these  measurements  were  used 
for  preliminary  design  considerations  and  are  compared  to  the 
orbits  which  were  obtained  by  numerical  integration  of  the  tra- 
jectory equation. 
1.   MOTION  OF  A  CHARGED  PARTICLE  IN  A  UNIFORM  MAGNETIC  FIELD 

A  longitudinal  type  homogeneous  magnetic  field  spectrometer 

has  commonly  been  used  in  the  analysis  of  p  spectra  in  nuclear 

U   5 
physics  '.   The  theory  of  motion  of  a  charged  particle  in  a 


15 


uniform   magnetic   field  was    originally  published  by    Busch      in    1926, 
upon   which   the  technique   was    based. 

The   trajectory   of   a    charged  particle    in    a    uniform   field    is    a 
helix,    where   the  particle,    after    one  complete   revolution    of  period 
T,    will    focus    again   on    the   axis    from  which    it    originated.      As    seen 
in   Figure   2,    the   distance  Z    ,    at   which   the    ion   returns    to   the  Z 
axis,    depends    on   the  momentum,    p,    the   angle,    0,    and   the   magnetic 
field,    B.      The   relationship   between   these   quantities    can   be   derived 
as    follows: 

Let   p   =   vector    momentum   of   an    ion   at    the  origin 

of   the    coordinate   system  at    t    =   0 
and   q   =    charge    of    ion c 
The  magnetic   field    is    described  by 

-♦AAA 
B   =    B      r    +    B  ,    e*    +   B      z 
r  0  z 

where        B   =  B,  =  0,  B   =  constant. 
x  p  3   z 

The  boundary  conditions  seen  from  Figure  2  are* 

at  t  =  0;  r  =  0,  z  =  0,  p   =  p  cos  0,  p   =  p  sin  0, 

and  at  t  =  T:  r  =  0,  z  =  z  . 

o 

In  order  to  get  the  relation  Z  =  f(p,0,B)  we  must  integrate 
the  force  equation,  F  =  q  v  x  B,  and  apply  the  boundary  conditions 
at  t  =  0,  and  t  =  T.   Thus 

F  =  q  v  x  B, 
and  in  component  form 

*  ■-  £l-r(if)2]  =«vA  =  WiH   (1) 
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F  ,  =  m 


•*)$■  -m 


q  v  B   = 
^   r  z 


{%)%    <2» 


F   = 


d  z  _  n 

m  — -  =  0. 

dt 


(3) 


Integrating  (3)  and  applying  the  boundary  conditions  at  t  =  0 

gives 

dz  „ 

m  -r—   =  constant  =  p  cos  Q 
dt  r 


z  =  |(p/m)  cos  0   t . 


(4) 


From  (2),  upon  multiplying  by  r  and  integrating 


m 


™m>  r2y 


^)Bz 


m 


dt 


r2/d^ 
(dt 


"qBz     d     ,    2, 

*      -2"   dt    (t     ) 


dt  2m 


=    constant    =   k . 


(5) 


Now  using    (5)    in    (1)    gives 


d2r  ,2  ,  n      .  2 

m -  mrk      =    -qrk   B      =    2rmk 

dt2 


2 

d  r  ,2 

— -  +   rk      =   0 

dt 


which  has  the  solution 


ikt    „    -ikt 
r  =  Cl  e     +  C2  e 


From    the    initial    conditions    at    t    =   0 


C      =    -C      =  p   sin    Q/2ikm 
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giving 


p  sin  0   .   ,  ^ 
■  ^ sin  kt«  (6) 


Now  eliminating  t  from  equations  (4)  and  (6)  gives 


p   s  in   0       .        .         km 

=  *H^; —  Sln  (: 


km  vp    cos    0 

and  at    t    =  T,    this    reduces    to 

p    cos    0 

Z      B      =    2tt  = .  (7) 

o      z  q  v     ' 

Therefore,  a  detector  at  Z   with  a  field  B  will  "see"  a  particle 

o  z              r 
— ► 

of  charge  q  scattered  with  momentum  p  at  an  angle  0  with  the  z 

axis.   Notice  that  there  are  two  ways  to  observe  different 

p  cos  0  values:  (1)  hold  B   constant  and  change  Z   by  moving  the 

detector,  (2)  fix  the  detector  at  a  particular  value  of  Z   and 

o 

vary  B  .   Thus  a  knowledge  of  equation  (7)  allows  one  to  position 

a  detector  to  collect  all  ions  within  a  given  momentum  interval 

I Ap|  at  |p|  scattered  into  a  conical  shell  AO  at  0.   In  order  to 

make  observations,  then,  one  would  simply  choose  the  values  of 

p  cos  0  to  be  observed,  solve  the  equation  for  various  values  of 

Z  and  B   that  will  allow  them  to  be  seen,  and  set  the  detector 
o       z 

and  field  accordingly. 

The  solid  angle  in  a  spectrometer  having  asrimuthal  symmetry 
depends  on  0  and  p  as  follows:   recall  that  dQ  =  2n  sin  0  dO 
where  we  have  integrated  over  d0.   An  expression  for  dO  can  be 
obtained  by  differentiating  equation  (7)  as  follows: 

Z   =  —r-  \   P  cos  0 
O    ^qB. 

and    since    B      =    constant,    then 
z  ' 

qB 


dZ      =  =-^    (cos    0   dp    -   p   sin   0  do) 

O  2TT        V  f  r-  / 
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qB       dZ 
,n    _         dp                        z          o 
and  dG     - pr    -    77— : rr     . 

p    tan    0         2tt  p    sin    0 


However 


giving 


3"  m     ^ 

p   =    (2mE)  2  or    dp    =    (-^)      dE 


He-  ^B      dZ~ 

do  dE 


2E   tan    0         ~     ] % 


2n(2mE)  "  sin   0 
Using   this    the   solid   angle    is    expressed   as 

(  2b  K2 

,  ~  ,TT      COS      9,        __  I    4         Z        \  ,  0, 

Therefore,  as  seen  from  equation  (8)  the  solid  angle  depends  on 
the  angle  of  scatter,  0,  the  trajectory  intercept  spread,  dZ  , 
allowed  by  the  detector,  and  the  energy  of  the  scattered  ion,  E. 

As  an  example,  consider  the  elastic  collision  He   on  He 
where  0  =  30°,  E  =  300  eV,  dE  =  5  eV,  B   =  200  gauss,  and 
dZ  =  -  3mm.   Then  dQ  =  O.058  steradians,  a  value  typical  of  ' 
solid  angles  for  the  nonuniform  magnetic  field  apparatus  actually 
used. 

2.   ANALOGUE  MEASUREMENT  OF  ORBITS 

Prior  to  obtaining  any  numerical  solutions  for  the  orbit 
equation,  the  technique  of  simulating  the  orbit  of  a  charged 
particle  in  a  magnetic  field  using  a  current  carrying  flexible 
wire  was  employed  to  obtain  preliminary  design  information  on  the 

various  experimental  parameters, 

7 
The  experiment,  carried  out  by  Kelly  ,  xs  based  on  the  fact 

that  a  current  carrying  wire  consists  of  charged  particles  moving 

with  some  velocity  along  the  wire,  causing  the  wire  to  experience 

a  v  x  B  force  in  a  magnetic  field,  the  same  type  of  force  charged 

particles  encounter. 
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Consider    a    small    length,    ds ,    of   flexible  wire   carrying  a 
current,     i,    and   under    a    tension,    T,    in   a    magnetic    field,    B. 


p  =  particle   momentum 

q  =  particle    charge 

p  =  radius    of    curvature 

m  =  mass/length 

i  =  current 

T  =  tension 

t  =  unit    vector    along 

tangent 
— • 

11  =  unit    vector    along   normal 


Then    since    the   system    is    in    equilibrium    the   sum   of    the   forces    is 

zero,    giving 

"•        ~*  ds  -*  ~* 

ids(T    x    B)    -   T   —  "H     -   mg    ds    =    0 

from  which    the    curvature   of   the   wire    is    shown    to   be 


T1         i     ,""•        -*.         mg 

_    =    —     (T     x     B)      -    — - 
P         T    ^  ;  T 


(1) 


To  get  the  path  of  a  charged  particle  in  a  magnetic  field,  recall 

that 


dp_  B 


dt 


q(v  x  B) 


~*    P    v 

but      t  =  *-  =  —  and  ds  =  vdt ,  for  constant  velocity,  v. 
p    v 

— »        — ► 

Therefore  v  -f~  =    vp   -3—   =    qv(T    x    B)  , 

ds  ^    ds         ^    v  ' 


however 


dT     _  T\_ 
ds    ""  P 


and  the  curvature  of  path  of  a  charged  particle  is 

J  =  a  (r  x  B)  . 

p   p  v 


(2) 
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Therefore   the  wire   and  particle  paths   will    coincide   when 

— * 
i    (T    x    B)     -    ^   =   £    (t    x    B).  (3) 

Now  for  our  case  the  magnetic  force  is  usually  100  times  the 
gravitational  force,  hence  the  gravitational  term  can  be  neg- 
lected and  equation  (3)  reduces  to 

?  =  £.  (4) 

Thus  when  equation  (4)  is  satisfied  it  allows  the  simulation  of 
charged  particle  trajectories  by  a  current  carrying  wire  under 
tens  ion . 

The  experimental  arrangement  used  to  make  the  wire  orbit 
measurements  is  shown  in  Figure  3.   An  uninsulated  wire  con- 
sisting of  l4  strands  of  number  44  copper  wire  weighing  0.24 
grams /meter  was  used,   Currents  of  up  to  4*5  amps  could  be 
passed  through  the  wire,  however  at  the  higher  currents  the 
wire  oxidized  and  became  less  flexible.   The  tension  in  the 
wire  was  normally  1000  dynes  and  was  measured  by  a  modified 
ammeter  movement  which  was  calibrated  against  a  set  of  weights. 
This  method  was  chosen  over  the  usual  pulley  arrangement  since 
the  frictional  effects  in  the  pulley  could  not  be  sufficiently 
reduced.   The  departure  angle,  0,  was  measured  to  an  accuracy  of 
+  30  minutes  of  arc  by  reading  the  longitudal  position  of  a  copper 
sleeve  at  the  instant  it  made  electrical  contact  with  the  wire. 
The  measured  departure  angle  was  then  corrected  to  account  for  the 
curvature  of  the  wire  between  the  teflon  ball  and  the  copper 
sleeve. 
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For    a    given  magnetic  field   and   distance   Z      the   current 

o 

through  the  wire  was  gradually  increased  until  the  desired  de- 
parture angle  was  obtained.   Then  the  tension  and  the  current  were 
measured  giving  the  momentum  (for  an  assumed  charge)  which  the  wire 
trajectory  represented.   The  shape  of  each  trajectory  was  obtained 
by  measuring  r  as  a  function  of  Z.   Thirty  seven  different  tra- 
jectories were  measured  and  tabulated  in  reference  7  which  also 
includes  a  detailed  discussion  of  the  method. 

A  comparison  of  these  measurements  with  the  computed  tra- 
jectories showed  that  the  measured  maximum  r  values  were  about 
17%  larger  than  the  calculated  maximum  r  values  and  the  measured 

Z  values  were  about  8%  larger  than  the  calculated  Z   values.   The 
o  o 

natural  stiffness  of  the  wire  is  believed  to  be  the  primary  cause 
of  the  discrepancy.   A  comparison  of  the  wire  orbit  and  the  calcu- 
lated trajectories  is  shown  in  Figure  4  for  an  H   ion  with  a  de- 
parture angle  of  48.17  degrees  and  scattered  energy  of  243  eV. 
3.   MAGNETIC  VECTOR  POTENTIAL 
The  trajectory  equation  has  as  one  of  its  parameters  the 
magnetic  vector  potential.   Therefore,  the  measured  field  of  the 

magnet  must  be  related  to  this  quantity  before  the  orbits  can  be 

— » 
found.   The  well  known  expansion  of  A  in  terms  of  the  field  is 

6 
found  by  following  the  development  used  originally  by  Busch  , 

Since  there  are  no  currents  in  the  region  of  interest 

V  x  B  =  0 
and,  consequently 

V  x  (V  x  A,  =0 
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is  the  equation  to  solve.   Now 

V  x  (V  x  A)  =  V(V  •  A)  -  V2^  (1) 

and  recalling  that  J   =   J^   &      only,  then  A  =  -, —  — * — *-  dv 

0  iiTT  J  |-»  I 

J   M 
requires  that  A  =  A,0  and  A   =  A   =0.   Also,  since  cylindrical 
symmetry  exists,  A,  is  not  a  function  of  0.      Thus  in  cylindrical 
coordinates,  (r,0,z),  equation  (;1)  becomes 

a_%  +  ii%^  +  aJi  =  o.  (2) 

-  2    r  Br      2   _  2  l  ; 

or  r     oz 

Now   assume   a    solution 


Art(r,z)    =  A(r,z)    =  E        a    (z)r    , 
0^     '     '  v     '     '         n=o      n 


then 


dA        v  n-l 

r —  -  I  na      r 
or  n 


-  =  Ln(n-l)a      r 

^2  v         '    n 

or 


o2A       r      "      n 
— —  =  E  a      r     . 
^    z  n 

oz 


Substituting    into   equation    (2)    gives 


v*       /       -.  n         n-2  ^  n-2     ^  n-2     ^      "      n 

E   n(n-l)a   r  +  E  na   r         -Ear         +  E  a      r      =0 

v         '    n  n  n  n 

E  [n(n-l)    +  n-1]  a    r    ~    +  E  a      r      =0 


E    (n-l)(n+l)    a    rn    2+  S  a"    r"    =   0, 

'  v  '       n  n 

Each  coefficient  must  be  zero  and  shifting  indices  in  the  second 

term  gives   »  oo 

E   (n+l)(n-l)a  tn~    +   E  a  jt    ""   =  0, 
n=0v    /v    y  n      n=2  n-2 
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and  for       n  =  0:   -l(aQ)  =  0;  a   =  0 

n  =  1:    0(a)  =  0;  a   =  arbitrary  function  of  z 

n  ^  2:    (n+l)(n-l)a  +  a   _  =  0, 
v    '  v    '  n    n-2 

The   recurrence   relation   for    the   coefficients    is 

it 

n-2 

a      = 


n         (n+l)(n-l) 

and   since   a„    =   0,    there   are   only    odd  n   terms    in   the   solution. 
The   general    term   is 

(-l)k  (2k) 

2k+1     '   22kk!(k+l)!      l 


This    gives 


»  „,    , ,  k  2k+l 


2(-D         ,*_.  (2k) 


A(r'z>    =S       kl   k+lH    (2>         ai 


.(k+1) 
k=0 


or 


3  5 

"    rJ        _(IV)    i 


"    r  , (IV)    r 

A(r,z)    =   axr    -    &1   -g-  +   a|       >  x_  +    ...     . 

Now   to   evaluate   the   coefficient   a    (z)    recall   from   B   =  V    x  A 

and   the   field    symmetry   that 

=    1  ii£Al  s   1    3     [a    ,2      «   r^  (IV)    r^   _  _ 

r      Br  r  or   L    1        al  X       al        192 


B 
z 


3  5 

,    »    rJ        ,    (IV)    rJ 

but    B    (r    =    0)    =    2a      giving   a^z)    =  -   B    (r    =   0). 


This    gives  2 

-  n  M  8       B>(r=0) 

Vr'Z>    ===0n}(nil)!    (2>  "^ 


(3) 
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Therefore,  by  knowing  B  (Z)  _0  we  also  know  A(-(r,z),  the  quantity 
used  to  solve  the  orbit  equation.   The  actual  measurements  were 
fitted  to  a  12    order  polynomial    in  z  such  that 

13    n  i 
Bz(z)  =  E  c/"1.  (4) 

With  this  expression  for  B  (z),  the  field  values  were  repro- 
ducible to  within  0.1%  of  the  experimentally  measured  values. 
Now  using  these  expressions  for  Acv(r,z)  and  B  (z)  in  the 

computer  program  one  can  calculate  the  trajectory  of  a  particular 

— * 

p/q  value  for  a  selected  magnetic  field. 

4.   ORBIT  EQUATION  IN  A  NONUNIFORM  FIELD 
The  thin  lens  spectrometer  with  its  axially  symmetric  non- 
uniform field  was  chosen  for  this  apparatus.   It  more  easily 
produces  the  large  field  required  to  bend  ions  and  it  gives  more 
spread  in  the  intercept  distances,  allowing  greater  resolution. 
The  equation  of  motion  for  a  charged  particle  in  an  axially  sym- 
metric nonuniform  magnetic  field  and  the  subsequent  employment 

of  this  field  geometry  as  a  p-ray  spectrometer  has  been  discussed 

8  Q 
by  several  authors  '  .   Our  derivation  of  the  trajectory  equation 

essentially  follows  that  of  Busch  ,  and  was  first  applied  to  our 

10 
particular  problem  by  Gagliano 

An  axially  symmetric  magnetic  field  can  be  represented  by 

the  curl  of  the  magnetic  vector  potential  having  only  an  azimuthal 

component,  A^  (see  previous  section).   The  motion  of  a  charged 

particle  in  such  a  field  is  governed  by  the  equation 

F  =  q  v  x  B  =  qvx(VxA). 
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In    cylindrical    coordinates,    (r,    0,    z),    the   force   components    are 

rd2r  ,d02  d0   ^rV 


dt 


dt 


dt      dr 


(1) 


F0   =  m 


fdr\/d0 
dt      dt 


+r 


d20 
dt2 


=  -q 


1    dr   ^rA0) 
r    dt  '    3r 


(2) 


d2z  d0   BA0 

Fz   =  m  ~2  =   *  r    dt  T£   ' 
dt 


(3) 


Now  multiplying   equation    (2)    by   r    and  using   the  relation 


d      _   dz    d  dr     S 

dt        dt  dz    +   dt   bx 


(4) 


gives 


dt   [mr      dt1    " 


-    q[ 


d(*A0) 
dt 


or 


Therefore 


d      r       2   d0  .      - 

dt  [rar      dT  +  q  rA0  ]    =   ° 


mr      -—  +   q  rArt   =    constant    of   the  motion, 
dt         ^        0 


or 


d0 
mr   -rr  +   q  Aw   =    constant/r  . 
dt  0 


At    r    =   0   the   left   hand   side    is    zero  requiring   that    the   constant 


be   zero.      Thus 


d0  / 

—  =    -    q  A0/mr . 


(5) 
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Now  substituting    (5)    into   equations    (1)    and    (3)    gives 


d  dr. 

dt    (m   dt> 


2  8Art 

s_  A    __1 

m       0     dr 


d     .      dz. 
-r—   (m  -tt)    - 
dt    v      dt ' 


2  SArt 

3_  A     _1 

m       0      dz 


Upon   rearranging    the  derivatives    in   these   two  equations    one   obtains 
2 


0  /  dz .       d   r         dr    r  dz     d     .      dz ,  -, 

m   "0  "^7  "  m(dT)    ~2  +  dz"  LdT  dl  (m  dT}  ] 

dz 


d  .      dr  - q 

dt v       dt '  *"        ™      ^v  ^+ ' 


and 


dz 


cSA 


— (m  — )    = 


0    _    dz  r  d  dz 


dt"       dt 


m       0      dz         dt      dz         dt 


r£<»27)  i- 


They    can   now   be    combined    into   the   one   equation 

2 


. dz .       d  r 
dz 


3  Art  2  9  Aw 

dz   L  m    a0     dz     '  m       0     Br 


=    0. 


(6) 


dz 
The   term    (-rr)      is    eliminated  using   the   vector    velocity   relationship, 

equation    (4),    and   equation    (5),    as    follows: 


-*>        dr    a  d0   *,        dz    a 

v   =  ■ —  r    +  r    ~*-  0   +  — —  z 

dt  dt    F        dt 


or 


2  2  2 

.dz,  2  .dr/         2    .d0. 


and   substituting   gives 


(£>2  =  [v2 -(>/]/[  1  +  (i)2i . 


dz 
Now  replacing    (-rr)       in   equation    (6)    leads    to  the   desired  trajectory 


equation 


dz2 


^+<£ )  i 


2         2      2 
where  K   =   p   /q    . 


[k2-a^-H 
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dA 


A 


0 


0     dz 


+    A0_£   =    O.       (7) 


Therefore  one  must  choose  the  value  of  p/q  to  be  used  and  know 
the  relation  for  A  =  f(B  ,z,r)  developed  in  the  previous  section 
in  order  to  integrate  the  trajectory  equation. 

The  trajectories  were  calculated  by  numerical  methods  on  the 

available  computer.   The  integration  routine,  originally  written 

10 
by  Gagliano   ,  used  a  fourth  order  Runge-Kutta  method.   However, 

it  was  rewritten  for  the  IBM  36O-67  computer  using  the  Hamming's 
modified  predictor -corrector  method  which  is  the  DHPGC  subroutine 
of  the  IBM  System/360  Scientific  Subroutine  Package.   The  com- 
plete computer  program  is  listed  in  Appendix  IV.   Figure  5  shows 
a  set  of  trajectories  calculated  on  the  computer  for   Li   +   Ne 
where  the  energy  of  the  incident  Li  ion  was  100  eV.   As  in  the 
case  of  the  uniform  magnetic  field  there  are  two  methods  of  de- 
tecting particles  scattered  at  different  angles:  (1)  by  moving 
the  detector  and  keeping  the  focusing  magnet  current  at  a  specific 
value;  (2)  by  keeping  the  distance  between  the  scattering  cell  and 
the  detector  fixed  while  varying  the  current  through  the  focusing 
magnet  coil.   Both  methods  were  used  and  acceptable  results  were 
obtained  from  each  method. 

5.   ANGULAR  AND  ENERGY  RESOLUTION 
As  is  shown  in  the  diagram  (Figure  6),  our  experimental 
arrangement  for  measuring  differential  scattering  cross  sections, 
da  (0),  as  a  function  of  scattering  angle  allows  a  spread  of 
angles  AG  at  0  to  be  admitted  by  the  detector.   Therefore  there 
is  some  angular  spread  associated  with  every  measurement  which 
causes  an  uncertainty  in  the  value  of  da  (0)  obtained  from  the 
measurement.   It  is  seen  from  the  sketch  (Figure  6)  that  the 

angular  spread  depends  on  the  size  of  AZ  ,  whose  value  is  determined 
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by  the  detector  apertures  and  9.   Thus,  AO  can  be  calculated  from 


the  relation  between  AO  and  AZ  .   This  relation  is  obtained  from 

o 


=  f (E,0,B,R' ,Z',C) 


(1) 


where 


E 

9 

B 

R< 

Z* 

C 


energy  of  the  incident  ion 

angle  of  scatter 

Magnetic  field  =  f(I);I  =  coil  current 

target  displacement  from  system  axis 

target  displacement  from  system  center 

a  constant  for  each  reaction  including  mass 

ratios,  inelastic  energy  losses,  etc. 


The  coil  current  is  constant  during  a  given  measurement  "making 
dl  =  0,  therefore 

';f 


:'5T' 


Also,    for    a   given    reaction,    C    is    a    constant   which      makes    dC   =   0. 
Thus    differentiating    (1)    leads    to 


az 


dz      = 
o 


dE+UZ°' 


*E/oR>z>~"Va9 


d0+ 


ER'Z 


dRf  + 


EOZ 


az 

r)Z 


dZ 


rE9R 


Solving    for    the   angular    spread  we  get 

^Z 


dO   =i 


J>0 
^Z 


dZ 


SO 


o      dZ 


'ER?Z 


?  -7  ? 


'ER*Z 


»7» 


dE+ 


ER'Z 


171 


dRc  + 


dZ 


EOZ1 


V 


3z 


dZ1 


(2) 


E^R 


The   values    of   the  various    quantities    in   the   expression    for 
the  angular    resolution  must    be    calculated  for    each   and  every    cross 
section    that    is   measured.      However    as    shown   by   example    in   Appendix 
I    the    terms    arising   from   the   energy    spread  and   the   finite    target 
size    increase   dO    by    less    than    10%   for    ion    beam   energy   spreads 
less    than    0„5   eV.      Therefore    the   angular    resolution    is 


*  -W       <*o 


(3) 
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when   the   terms    containing   dE,    dR*,    and   dzr    are    ignored. 

The    second    quantity   of    interest    is    the   spread    in   energies, 
dE    ,    that   particles    scattered   at   0  may   have   and   still   be   detected, 
As    can   be   seen   from  Figure   6,    this    quantity   depends    upon    AO   at    0 
and   also  upon    the   energy   of    the    incident    ion,    E,    and   the  reaction 
parameter    C.      Therefore 

E      =    E       (E,'0,C) 
s  s 

and  upon    differentiating 

/dE  \  (bE 

If   one   considers    only    elastic  processes,    where   the    incident    ion 
of  mass    m   interacts    with  a    target    of  mass   M,    then,    due    to   the 
conser vat  ion   of   energy   and  momentum      , 


2  2  2  2  2222 

„  „   M      -   m      +   2m      cos    0   +    2m   cos    0    \M      -   m      +   m      cos    0 

E      =    E 


(M    +  m)2 


Using    this 


and 


r3E  \ 

pr)0 '-  ^ 


dE  \  -2m  E      Sin   0 

s  s 


S0    /  (— 


-i 


-   m      +   m      cos    0 


giving 


d0    .  (5) 
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Using  equations  (3)  and  (5)  one  can  evaluate  the  angular  resolution 

for  a  particular  cross  section.   As  an  example,  for  the  reaction 

6Li+  on  He  at  300  eV,  I  =  13-0  amps,  dE  =  0,25  eV  and  0  =  39.6°, 

one  obtains 

dQ  =  o.38° 


and 


dE   =  3.77  eV 
s 


which  are  representative  values  of  the  reactions  that  were 
observed. 

6.   SOLID  ANGLE 
The  major  advantage  of  this  method  of  measuring  differential 
scattering  cross  sections  is  the  greatly  enhanced  solid  angle 
obtained  by  collecting  all  particles  scattered  into  a  conical 
shell  dQ  at  0.   The  solid  angle  for  this  system  is 

dQ  =  2tt  sin  0   dO . 

As  can  be  seen  from  Figure  6  all  the  particles  scattered  at  0 

must  cross  the  axis  in  AZ   at  Z   in  order  to  be  detected.   This 

o      o 

AZ  ,  then,  determines  d9  at  9  and,  consequently,  the  value  of 

the  solid  angle  at  0.   In  the  discussion  of  the  angular  resolution 

the  expression 

was  obtained.   This  is  the  relation  used  to  calculate  the  solid 
angle.   Thus 


dQ  (0)  a  2tt  sin  0  U 


SO 


Z 


dZ 
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However  only  290   of  the  available  36O   of  the  axial  angle  were 

open  because  of  the  way  the  scattering  cell  is  constructed  (see 

Chapter  III,  Section  5,  for  a  full  explanation).   Correcting  for 

this  and  expressing  d©  /Sz   in  degrees/centimeters,  the  solid  angle 

o    o 

expression  is  written 


|Q(0)=  O.O873  sin  0(||-j   dZQ 


d0(0)  =  0.0873  sin  ©  (  s-z-\      dz_  steradians. 

E 
This  expression  must  of  course  be  evaluated  for  every  trajectory 

to  be  measured  since  it  depends  not  only  on  the  angle  of  scatter, 

0,  but  also  on  AZ  and  the  derivative  (d©/dZ  )„.      All  of  these 
'  o  v      o'E 

quantities  depend  upon  the  particular  ion  trajectory  considered. 
Some  representative  values  of  the  solid  angle  as  a  function 
of  0  at  various  focusing  magnet  currents  are  shown  in  Figure  7. 
These  curves  are  for  the  reaction   Li   +  "  Ne  at  100  eV  incident 
ion  beam  energy.   Notice  the  dependence  of  the  solid  angle  on  the 
angle  of  scatter  and  the  focusing  magnet  current  (magnetic  field) . 
When  the  magnet  current  is  held  constant  the  intercept  distance 
changes  for  each  angle  of  scatter  causing  the  solid  angle  to  vary. 
These  values  at  the  solid  angle  are  typical  for  all  reactions 
observed. 
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Figure  7.   Solid  Angle  vs  Angle  of  Scatter  for  Elastic  Scattering 
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of  Li   on  Ne  at  100  eV  Incident  Ion  Energy 


38 


CHAPTER  III 
EXPERIMENTAL  APPARATUS 
1,   VACUUM  SYSTEM 

The  vacuum  shell  shown  in  Figure  8  consists  of  aluminum 
cans  with  neoprene  O-rings  at  the  connecting  flanges,  feed 
throughs,  and  end  plates  * 

Two  six-inch  diffusion  pumps  are  used  to  maintain  the 
vacuum.   A  National  Research  Corporation  Model  VHS6  diffusion 
pump  having  a  2400  liter  per  second  pumping  speed  is  connected 
to  the  detector  end  of  the  system.   Water  cooled  baffles  and  a 
liquid  nitrogen  cold  trap  are  located  between  the  detector  chamber 
and  the  diffusion  pump c   The  forepump  is  a  CENCO  HYVAC  45  mechan- 
ical pump  with  a  280  liter  per  minute  capacity.   A  stacked  half- 
moon  baffle  is  placed  in  the  foreline  to  reduce  backstr earning 
from  the  forepump.   At  the  opposite  end  of  the  system,  where  the 
ion  source  and  the  mass  analyser  are  located,  a  Dresser  Model 
DPD-6  diffusion  pump  is  utilized.   The  Dresser  pump  has  an  1800 
liter  per  second  pumping  speed  and  is  connected  to  a  Welch  Duo- 
Seal  forepump*   Water  cooled  baffles  and  liquid  nitrogen  trap  are 
installed  between  the  diffusion  pump  and  the  mass  analyser  chamber* 

Dow  Corning  704  diffusion  pump  oil  is  used  in  both  diffusion  pumps. 

-7 
This  arrangement  provides  ultimate  pressures  of  2.0  x  10    Torr 

-7 
in  the  detector  chamber,  and  1.0  x  10    Torr  in  the  ion  source 

chamber .   With  gas  in  the  scattering  cell  at  a  pressure  of 

-3 
2»0  x  10   Torr  the  pumps  maintain  the  detector  chamber  at 

2c0  x  10  "  Torr  and  the  mass  analyser  chamber  at  4.0  x  10    Torr. 

These  pressures  are  low  enough  to  insure  that  reactions  outside 
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Figure  8.   Experimental  Apparatus 


4o 


the  scattering  cell  are  negligible.  Tests  also  showed  that  the 
pressure  in  the  detector  chamber  is  low  enough  to  eliminate  any 
pressure  effects  on  the  electron  multiplier  which  served  as  the 
detector  for  the  scattered  particles. 

2.   ION  SOURCES 
In  order  to  generate  the  ion  beams  used  in  the  differential 
scattering  cross  section  measurements}  two  different  types  of  ion 
sources  were  used:  Lithium  ions  were  produced  by  a  thermal  emitter; 
helium  and  hydrogen  ions  by  an  electron  impact  source.   The  per- 
tinent facts  for  both  types  of  sources  are  reported  here  while  a 
more  detailed  account  of  our  developmental  work  on  the  electron 
impact  source  is  presented  in  Appendix  VI. 
Lithium  Ion  Source 

The  lithium  ion  source  used  was  essentially  the  source  de- 

12 
scribed  by  Heinz  and  Reaves   and  is  shown  in  Figures  9  and  10. 

It  is  a  thermal  emitter  consisting  of  a  very  porous  tungsten 
matrix  into  which  the  minerals  beta-eucrypt ite  (Lio0»Al_0  *2SiO  ) 
or  spodumene  (Li  O'Al  0  "4si0  )  have  been  fused.   The  lithium 
ions  are  produced  when  the  tungsten  matrix  is  heated  to  tem- 
peratures in  excess  of  1100  C  by  a  separate  filament  (see  Figure 
10) ,   The  emitter  was  enclosed  in  a  heat  shield  to  keep  filament 
power  requirements  low,   The  shielded  emitter  was  mounted  on  a 
8-1/2  inch  diameter  stainless  steel  flange  using  a  stainless  steel 
mounting.   In  front  of  the  emitter  was  an  extraction  electrode 
and  a  three  electrode  einzel  lens.   To  extract  the  ion  beam,  the 
emitter,  shield,  and  filament  were  placed  at  a  potential  corre- 
sponding to  the  desired  beam  energy.   The  filament  power  was 
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Figure   9.    Schematic    of   the  Lithium   Ion   Source. 
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o 
adjusted  for  an  emitter  temperature  of  about  1170  C  which  was 

found  to  be  an  optimum  operating  point.   The  effect  of  the  emitter 
temperature  on  the  total  ion  emission  is  shown  in  Figure  11,   The 
extraction  grid  was  grounded  producing  a  beam  of  the  desired  energy. 
The  einzel  lens  apertures  were  1,4  centimeters  in  diameter  and  the 
lenses  were  1,4  centimeters  apart.   The  center  lens  was  4.0  centi- 
meters from  the  extraction  electrode.   These  dimensions  were  ob- 
tained from  the  einzel  lens  relations  given  in  Appendix  VI  and 

-7 
produced  ion  beams  up  to  1.0  x  10    amps  in  the  scattering  cell 

when  the  emitter  was  operated  at  optimum  output.   The  ion  beam  ob- 
tained from  the  source  was  found  to  have  an  energy  spread  of  no 
more  than  0.5  eV  in  the  energy  range  15  to  800  eV , 

This  type  of  ion  source  has  several  advantages  over  plasma 
type  sources:  simplicity  of  operation,  small  gas  load  on  the  system, 
stable  and  highly  monoenergetic  beams  of  high  purity.   The  major 
drawback  of  course  is  a  limited  number  of  species  which  can  be 
produced  by  this  technique. 
Electron  Impact  Ion  Source 

This  ion  source  is  one  of  a  large  family  of  sources  which 
produces  ion  beams  by  creating  a  plasma  in  the  appropriate  gas 
and  extracting  ions  from  the  plasma  by  externally  applied  elec- 
trostatic f ields o   In  the  particular  source  shown  in  Figure  12, 
the  plasma  is  produced  by  electron  bombardment  of  the  low  pres- 
sure gas.   Electrons,  emitted  from  the  filament,  accelerated 
radially  outward  through  the  concentric  grid  and  decelerated  as 
they  approach  the  cylindrical  walls,  oscillate  between  the  fila- 
ment and  the  wall  until  they  are  collected  on  the  grid.   The 
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Figure  12.  Electron  Impact  Ion  Source, 


plasma  forms  between  the  grid  and  filament  and  the  ions  are  ex- 
tracted axially  from  it.   To  produce  an  ion  beam,  source  gas  is 

introduced  into  the  source  can  at  a  rate  that  will  maintain  a 

-3 
pressure  of  about  1  x  10   Torr  within  the  source.   The  potentials 

are  set  at  the  desired  ion  beam  energy  and  the  filament  power  is 
increased  until  enough  electrons  (usually  about  90  ma)  are  pro- 
duced to  maintain  the  discharge.   The  grid  is  usually  100  volts 
above  the  common  can  and  filament  potentials.   An  einzel  lens 

placed  directly  in  front  of  the  extraction  hole  helps  produce 

-8 
ion  beams  up  to  3.°  x  10   amps  in  the  scattering  cell.   Our 

measurements  indicate  a  rather  large  energy  spread  of  about  8  eV 

although  other  workers  using  the  same  source  have  reported  energy 

spread  as  low  as  .5  eV  . 

The  source  is  fairly  easy  to  operate  and  maintain,  although 
requiring  more  maintainance  than  the  lithium  source.   The  primary 
drawbacks  are  the  gas  load  it  imposes  on  the  vacuum  system  and  the 
fairly  large  energy  spread. 

A  more  detailed  account  of  our  developmental  work  on  this 
ion  source  is  presented  in  Appendix  VI. 

3.   MASS  ANALYSER 

The  ion  sources  used  tend  to  produce  a  variety  of  ions  from 
the  same  source.   This  was  especially  true  when  using  a  plasma 
type  ion  source.   Although  the  source  parameters  could  be  varied 
to  enhance  the  production  of  one  particular  ion,  other  ions  could 

^Private  communication  from  Dr.  W.  Aberth,  SRI. 
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not  be  reduced  to  insignificant  values.  Therefore,  it  was  decided 
that  mass  analysis  of  the  ion  beam  prior  to  its  entering  the  scat- 
tering cell  was  necessary. 

The  mass  analyser  shown  in  Figure  13  was  designed  and  cali- 

13  o 

brated  by  Strohsahl   .   It  has  a  45   angle  of  bend  in  order  to 

allow  for  both  high  beam  transmission,  up  to  90%?  and  a  resolution 

Am/m  of  about  40  allowing  mass  analysis  of  singly  charged  particles 

up  to  mass  number  40  and  energies  up  to  one  kilovolt.   Using  a  coil 

current  of  8  amps,  it  had  a  maximum  field  strength  of  about  4000 

/  o 

gauss  confxned  between  40   wedge  sector  shaped  pole  pxeces  with  a 

3/4  inch  gap  width.   The  pole  pieces  were  external  to  the  vacuum 
system  which  allowed  increased  flexibility  and  enhanced  the  vacuum 
integrity  of  the  entire  system. 

One  effect  resulting  from  the  use  of  the  mass  analyser  was 
the  reduction  of  background  signal  at  the  detector  by  a  factor 
of  100  or  more.   This  fact  made  it  worth  while  to  use  the  mass 
analyser  even  in  cases  where  the  ion  source  did  produce  only  one 
type  of  ion. 

4.   FOCUSING  MAGNET 

The  magnet  cell  consists  of  five  concentric  brass  spools 

wound  with  a  total  of  793&  turns  of  number  15  wire  connected  in 

such  a  way  (see  Table  I)  that  the  total  resistance  was  approximately 

10,75  ohms.   A  single  layer  of  1/4  inch  outside  diameter  copper 

tubing  was  placed  over  the  outermost  layer  of  wire  on  each  spool. 

Water  was  run  through  the  tubing  for  cooling,   Thus  it  was  possible 

to  use  up  to  14  amps  through  the  coils.   A  Sorenson  Nobatron  model 

DCR  150-15  power  supply  was  used  to  supply  the  required  magnet  coil 

current  regulated  to  +  0.1%,. 

48 


Ion  Source  Flange 

Thermal    Cement 
Cooling  Coils 

Ion  Source  Spool 
Vacuum  Weld 


Vacuum  Weld 


Lucite   Window 


Pumping  Flange 
Adapting  Flange 


0        1        2        3       4       5* 

Ti  1 1 1       i       i       i       i       i 
SCALE   INCHES 


Figure    13.    Schematic   of  Mass    Spectrometer. 
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TABLE  I 

FOCUSING  MAGNET  DATA 

i 
ITEM              | 

1 

Number  of  layers  of    : 
winding                j 

'coil 

1 

COIL    COIL 

COIL 

4 

COIL   j TOTAL 

__5  _  _J_ 

1. 

22 

18 

16 

12 

14 

82 

2 

Average  number  of 
turns  per  layer*       1 

97 

95 

97 

99 

98 

3. 

1 
Total  number  of        \ 

turns  on  spool 

(item  1  x  item  2) 

i 

2116 

1710 

1552 

1188 

1372 

7938 

4. 

! 

Measured  resistance 

i 
( ohms ) 

20.46 

18.36  20.86 

18.04 

23.92 

10.75 

:  5. 

Length  of  wire  (ft) 
(item  4  x  1000/2.58)**  j 

6010 

6790   7560 

6790 

9300 

36,450 

6. 

Inside  diameter  of 
spool  (inches) 

7.750 

11.155  14.300 

17.250 

19.750 

7.750 

7. 

Inside  diameter  of 
winding  (inches) 

8.000 

11.400  14.550 

17.500 

20.000 

8.000 

8. 

Outside  diameter  of 
spool  (inches) 

J 

11.140 

14.285 

17.235 

19.735 

22.250 

22.250 

*Based  on  actual  count  of  several  random  layers. 


** Based  on  nominal  resistance  for  15  wire  of  3-25  ohms/1000  ft 
at  25°C. 
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Since  the  magnetic  vector  potential  can  be  expressed  in  terms 
of  the  axial  component  of  the  magnetic  field  on  the  axis  only  this 

quantity  was  carefully  measured.   The  measurements  were  carried 

7  .       . 

out  by  Kelly   using  a  Hall  probe  calibrated  in  a  magnetic  field 

varied  from  0  to  5  kilogauss  against  a  known  nuclear  magnetic 

resonance  probe  accurate  to  +  1  gauss  at  1000  gauss.   Figure  14 

illustrates  the  method  used  to  position  the  magnetic  field  probe. 
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Figure  14.  Arrangement  for  Measuring  the  Magnetic  Field. 

The  probe  was  placed  inside  an  aluminum  pipe  on  the  tip  of  a 
wooden  rule  accurate  to  0.5  mm.   The  wooden  rule  provided  an 
adequate  measurement  of  the  distance  of  the  Hall  probe  from  the 
center  of  the  magnet.   The  aluminum  pipe  holding  the  probe  was 
fitted  through  templates  at  either  end  of  a  larger  aluminum  pipe 
held  rigid  inside  the  magnet  coil  and  positioned  concentric  with 
the  magnet  coil  axis.   The  end  plates  could  be  rotated  to  different 
radial  angles  and  they  had  seven  holes,  1.27  cm  apart,  along  a 
radius  which  allowed  field  measurements  at  various  r  and  0  values. 
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The  axial  component  of  the  magnetic  field  on  the  axis  was 
measured  at  4  centimeter  intervals  at  current  settings  of  294,6, 
8,10,12,  and  l4  amps.   A  plot  of  these  measurements  is  shown  in 
Figure  15«   These  values  were  fitted  to  a  12 —  order  polynomial 
in  Z  and  used  to  calculate  the  magnetic  vector  potential  as  ex- 
plained in  Section  3  of  Chapter  II.   Figure  16  shows  the  value 
of  the  axial  component  of  the  magnetic  field  measured  at  various 
r  and  0  values  at  the  center  of  the  coil.   These  and  other  mea- 
surements tabulated  in  reference  7  demonstrated  that  the  magnetic 
field  was  symmetric  about  both  the  center  of  the  coil  and  the 
coil  axis  to  within  0.1%.   This  allowed  use  of  the  coil  axis 
and  coil  centerline  as  references  for  alignment  of  the  entire 
system  to  the  magnetic  field  axis. 

5.   SCATTERING  CELL 

The  scattering  cell,  shown  in  Figures  17  and  18,  constructed 
of  stainless  steel,  is  designed  to  allow  particles  scattered  at 
all  a^imuthal  angles  to  leave  the  scatter  ing  cell  and  eventually 
be  collected  at  the  detector „   Figure  18  illustrates  how  this  was 
accomplished  by  opening  the  scattered  beam  exit  slit  about  290 
around  the  scattering  cell. 

Attached  to  the  front  of  the  scattering  cell  is  a  circular 
aperture  of  variable  diameter  which  can  be  adjusted  from  outside 
the  vacuum  system.   Through  appropriate  wiring  the  ion  currents 
to  the  front  section  of  the  scattering  cell  can  be  monitored,   In 
the  front  section  of  the  scattering  cell  a  0.5  centimeter  diameter 
by  1.0  centimeter  long  beam  entrance  tube  was  drilled.   These 
dimensions  were  chosen  to  reduce  gas  flow  out  from  the  scattering 
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Figure  15.  Axial  Magnetic  Field  Strength  vs  Distance. 
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cell  and  allow  maximum  flexibility  in  the  1.0  centimeter  diameter 
by  1.0  centimeter  long  scattering  chamber  itself.   Into  the  bottom 
of  this  chamber  a  1/8  inch  outside  diameter  copper  tube  was  placed, 
through  which  target  gas  entered.   A  precision  gas  metering  valve 
placed  at  the  other  end  of  the  inlet  tube  was  used  to  allow  pre- 
cise control  of  the  gas  flow  into  the  scattering  chamber.   A 
VG-1A  ion  gauge  was  located  so  that  its  entrance  tube  was  placed 
90   from  the  gas  inlet  tube.   The  narrowest  portion  of  the  path 
to  the  ion  gauge  was  0.5  centimeters  in  diameter,  occur ing  where 
it  entered  the  scattering  chamber  itself.   From  the  rear  of  the 
scattering  chamber  the  front  section  of  the  scattering  cell 
sloped  upward  to  the  outer  edge  of  the  cell  at  an  angle  of  49 
from  the  horizontal.   This  surface  acts  as  the  front  edge  of  the 
scattered  beam  exit  slit.   The  back  edge  of  the  slit  was  formed 
by  the  face  of  the  rear  of  the  scattering  cell  which  sloped  upward 
at  an  angle  of  36  .   As  shown  in  Appendix  III  these  two  angles 
determined  the  gas  target  thickness.   The  inside  of  the  rear  of 
the  scattering  cell  was  opened  up  to  a  2,54  centimeter  diameter 
in  order  to  allow  free  travel  of  the  unscattered  ion  beam  to  the 
beam  collector.   This  after  section  is  also  wired  so  that  ion 
currents  to  it  can  be  measured.   The  rear  section  of  the  scat- 
tering cell  can  be  moved  horizontally  from  outside  the  vacuum 
system.   In  this  way  the  target  thickness  could  be  varied.   The 
opening  could  be  adjusted  to  any  desired  length  from  fully  closed 
to  up  to  0.8  centimeters  when  fully  opened. 
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The  front  and  rear  sections  of  the  scattering  cell  were  placed 
inside  a  stainless  steel  tube  (0.l8  centimeter  wall  thickness)  and 
separated  from  it  by  teflon  insulation.   The  support  tube  ran  the 
length  of  the  entire  scattering  cell  but  was  cut  away  at  the  exit 
slits  so  that  only  19%  of  the  tube  was  left  to  obstruct  the  scat- 
tered beam.   This  area  was  accounted  for  in  the  calculation  of 
the  solid  angle.   The  tube  has  four  variable  length  legs  on  each 
end  (a  total  of  eight).   The  individual  legs  in  each  set  are 
placed  90   apart  around  the  circumference  of  the  tube  end.   This 
arrangement  gives  very  stable  support  to  the  cell  and  yet  allows 
easy  adjustment  for  alignment. 

The  beam  collector  consists  of  three  grids  and  a  collector. 
Normally,  the  grid  farthest  from  the  collector  was  grounded,  the 
middle  grid  was  given  s  small  positive  potential  to  keep  slow 
positive  ions  from  reaching  the  collector,  and  the  grid  nearest 
the  collector  was  given  a  small  negative  potential  to  suppress 
secondary  electron  emission  from  the  collector.   The  collector 
was  kept  at  ground  potential.   The  beam  collector  and  its  grids 
are  placed  inside  a  metal  cup  which  both  shields  them  from  ex- 
ternal charged  particles  and  also  acts  as  a  gas  cap  over  the  rear 
end  of  the  scattering  cell.   The  beam  collector  assembly  is 
mounted  on  a  rod  which  allows  the  assembly  to  be  pulled  away  from 
the  rear  of  the  scattering  cell  and  then  swung  out  of  the  path  of 
the  incident  ion  beam.   This  is  done  to  allow  alignment  of  the 
entire  apparatus  along  the  magnet  axis  and  to  allow  the  incident 
ion  beam  to  fall  directly  on  the  detector  for  calibration  purposes. 
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6.   SCATTERED  ION  DETECTOR 
The  differential  scattering  cross  sections  to  be  measured  are 
extremely  small,  and  therefore  it  will  be  necessary  to  use  some 
sort  of  amplification  so  that  the  scattered  particle  current  can 
be  measured  on  an  ammeter.   A  Model  306  Magnetic  Electron  Multi- 
plier manufactured  by  Bendix  Corporation,  Cincinnati,  Ohio,  was 

chosen  as  the  beam  detector  and  amplifier.   This  multiplier,  de- 

14 
veloped  by  Goodrich  and  Wiley   ,  could  provide  reproducible  cur - 

rent  gains  up  to  10  ,  sufficient  for  the  projected  measurements. 

The  multiplier,  shown  schematically  in  Figure  19,  consists 

of  two  plane  parallel  glass  surfaces  each  5>5   centimeters  long 

and  placed  0.5  centimeters  apart.   The  glass  strips  are  coated 

on  their  facing  sides  with  a  thin,  high  resistance  (approximately 

o 

10   ohms)  layer  of  tin  oxide  and  antimony.   Attached  to  the  top 
of  the  field  strip  is  an  entrance  grid  of  etched  nickel  mesh,  1.9 
centimeters  wide  and  1.2  centimeters  high,  with  greater  than  90% 
transparency.   Attached  to  the  top  of  the  dynode  strip  and 
directly  behind  the  entrance  grid  is  a  tungsten  cathode  2.54 
centimeters  long.   Permanent  Alnico  V  magnets  placed  outside  both 
the  field  strip  and  the  dynode  strip  produced  a  uniform  magnetic 
field  of  400  to  600  gauss  in  the  space  between  the  strips.   This 
magnetic  field  is  oriented  perpendicular  to  the  long  dimension  of 
the  strips  and  parallel  to  their  faces.   At  the  bottom  of  the 
strips  and  between  them  a  stainless  steel  anode  is  located.   The 
entire  structure  is  quite  rugged  and  easy  to  handle.   Its  outside 
dimensions  are  10.2  centimeters  high,  3»35  centimeters  long,  and 
2.06  centimeters  deep. 
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Figure  19.   Operational  Schematic  of  the  Bendix  Model  306 
Magnetic  Electron  Multiplier. 
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In  the  multiplier's  normal  mode  of  operation  a  potential 
difference  is  applied  the  length  of  each  strip  in  such  a  way  that 
the  end  of  the  strips  nearest  the  anode  is  the  more  positive  end, 
and  so  that  the  field  strip  is  positive  with  respect  to  the  dynode 
strip.   This  arrangement,  along  with  the  magnetic  field  of  the 
permanent  magnets,  produces  crossed  electric  and  magnetic  fields 
in  the  region  between  the  two  surfaces.   In  this  field  configu- 
ration the  motion  of  any  electron  released  from  rest  at  any  point 

— ►   — » 
on  the  dynode  strip  is  a  cycloidal  motion  in  an  E  x  B  direction 

(down  the  long  dimension  of  the  strip).   Because  the  equipotential 
lines  are  slanted  with  respect  to  the  dynode  strip  surface  the 
electron  strikes  the  surface  before  the  completion  of  the  first 
cycle.   For  sufficiently  high  applied  potentials  the  electron 
has  enough  energy  to  cause  secondary  emission  with  a  secondary 
emission  ratio  greater  than  unity.   The  secondary  electrons  thus 
generated  repeat  the  process  of  the  original  electron,  until  most 
all  electrons  reach  the  anode  and  are  collected.   Therefore  any 
current  incident  on  the  tungsten  cathode  is  read  out  on  the  am- 
meter connected  to  the  anode  multiplied  by  some  measurable  gain. 

As  the  scattered  ion  detector  in  the  experimental  apparatus, 
the  multiplier  has  to  be  operated  in  a  magnetic  field  of  up  to 
200  gauss.   The  effect  of  this  magnetic  field  causes  the  secondary 
electrons  of  the  multiplier  to  drift  perpendicular  to  the  long 
axis  of  the  dynode  strip.   Thus  some  of  these  electrons  drift  off 
the  side  of  the  dynode  strip  and  never  reach  the  anode,  thereby 
reducing  the  gain  of  the  multiplier.   The  rate  of  drift  increases 
with  increase  in  the  magnetic  field,  causing  the  multiplier  gain 
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to  decrease  with  increasing  field.   Shielding  the  multiplier  from 
the  external  field  was  attempted  using  p.  -metal,  netic,  and  conetic 
foils.   This  proved  unsatisfactory  for  two  reasons.   First,  the 
large  amount  of  shielding  material  required  so  disturbed  the 
focusing  field  that  scattered  particle  trajectories  were  not  pre- 
dictable.  Second,  the  gain  of  the  multiplier  could  not  be  re- 
produced whenever  the  focusing  field  was  changed  and  then  re- 
turned to  the  original  value.   This  was  probably  due  to  the  hy- 
steresis in  the  shielding  material.   The  most  satisfactory  system 
has  been  to  operate  the  multiplier  inside  an  aluminum  shield  to 
protect  it  from  random  background  and  to  choose  the  various  strip 
potentials  which  give  the  highest  gain  and  lowest  noise.   For  a 
particular  set  of  potentials  the  multiplier  gain  was  measured  at 
different  focusing  magnetic  field  values.   This  was  done  by  varying 
the  magnet  coil  current  and  keeping  the  multiplier  at  a  fixed 
distance  from  the  magnet  center  and  also  by  fixing  the  coil  current 
and  varying  the  distance  of  the  multiplier  from  the  magnet  center, 
both  methods  gave  essentially  the  same  gain  for  equal  values  of  the 
axial  magnetic  field.   A  plot  of  these  measurements  and  a  10 — 
order  fit  of  the  data  points  is  shown  in  Figure  20.   The  fitted 
value  was  used  in  the  actual  calculation  of  cross  sections.   An 
actual  gain  measurement  is  made  as  follows:  with  the  multiplier 
installed  inside  the  apparatus,  the  focusing  magnet  current  and  the 
detector  distance  are  established  at  their  desired  values.   Then 
the  field  strip  and  the  anode  of  the  multiplier  are  grounded  while 
the  cathode  is  connected  to  an  ammeter.   The  ion  beam  from  the 
source  is  directed  onto  the  cathode.   This  current  (which  has  to 
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be  kept  below  10   "  amps  to  prevent  saturation  of  the  multiplier 
output)  is  then  measured.   Then,  without  changing  the  current 
incident  on  the  multiplier  cathode  the  multiplier  is  placed  in 
its  normal  operating  condition  and  the  anode  current  is  determined. 
The  ratio  of  the  anode  current  to  the  ion  beam  current  measured 
at  the  cathode  is  the  gain.   This  process  is  then  repeated  for 
each  new  set  of  magnet  currents  and  detector  distances. 

Measurements  showed  that  the  gain  did  not  vary  with  gas 

-4 

pressure  up  to  2.0  x  10   Torr ,  a  pressure  much  higher  than  the 

multiplier  ever  encountered  when  data  was  being  taken.   The  gain 

.6    .7     4 
was  measured  using  Li  ,  Li  ,  He  ,  H  ,  and  H  ions  and  no  ion  or 

mass  dependence  was  detected.   The  gain  was  found  to  be  inde- 
pendent of  the  anode  current  up  to  10    amps.   Above  these  cur- 
rents the  multiplier  tends  to  saturate  since  the  anode  current 
is  approaching  the  strip  current  in  magnitude.   To  prevent  slow 
ions  from  entering  the  detector,  a  grid  was  placed  over  the  de- 
tector shield  aperture.   Normally,  a  small  positive  potential 
(about  20  volts)  was  sufficient  to  screen  out  this  background. 
The  detector  signal  falls  about  50%  as  the  grid  voltage  is  in- 
creased from  0  volts  to  about  10  volts.   As  the  voltage  is  further 
increased, the  detector  signal  remains  essentially  constant  until 
the  scattered  ion  beam  energy  is  reached,  then  the  detector  sig- 
nal again  decreases  sharply.   The  grid  also  permits  a  crude  anal- 
ysis of  the  scattered  ion  energy. 

After  about  two  months  operation  the  multiplier  strips  be- 
came coated  with  vacuum  pump  oil.   This  caused  large  increases  in 
the  multiplier  noise  resulting  in  false  current  measurements. 
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Cleaning  the  strips  with  a  pencil  eraser  to  remove  the  deposits 
reduced  the  noise  and  did  not  effect  the  multiplier  gain. 

Subject  to  the  limitations  outlined  above  the  magnetic  elec- 
tron multiplier  proved  to  be  a  reliable  detector  for  the  scat- 
tering apparatus.   Although  initial  measurements  of  the  gain  were 
tedious,  once  they  were  made,  subsequent  checking  showed  the  gain 
to  be  stable. 

7.   SYSTEM  ALIGNMENT 

In  order  to  detect  particles  scattered  at  a  particular  angle 
the  detector  must  be  placed  on  the  focusing  magnetic  field  axis 
at  the  point  where  the  scattered  particles  focus  on  the  axis. 
The  program  for  calculating  the  intercept  distance  assumes  that 
scattering  took  place  on  the  axis  at  the  center  of  the  focusing 
magnet  field.   Therefore  the  accuracy  of  the  angle  of  scatter 
determination  depends  on  how  well  the  scattering  center  and  the 
detector  are  aligned  to  the  magnetic  field. 

Since  the  geometric  center  line  of  the  focusing  magnet  coil 
was  found  to  be  the  magnetic  field  axis  to  within  0.1%,  this 
geometric  centerline  was  used  as  the  reference  for  aligning  all 
other  components  to  the  magnetic  field  axis.   To  perform  the 
alignment,  crosshairs  were  placed  on  either  side  of  the  magnet 
spool  to  mark  its  axis.   Then  a  small  laboratory  laser  was 
positioned  so  that  its  beam  lay  on  the  magnet  axis  established  by 
the  two  crosshairs.   The  crosshairs  were  removed  and  the  line 
established  by  the  laser  beam  was  used  as  the  reference  for 
aligning  all  other  components  to  the  focusing  magnet  axis.   This 
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method  made  possible  the  aligning  of  the  scattering  cell,  the 
detector,  and  the  focusing  lenses  to  the  magnetic  field  axis  to 
an  accuracy  of  1.0  milimeter  radical  displacement  off  the  axis. 
Analysis  of  computed  trajectories  showed  this  error  to  be 
essentially  undetectable. 

The  incident  ion  beam  was  aligned  by  adjusting  the  position 
of  the  source  and  the  mass  analyser  with  the  ion  beam  on  until  the 
maximum  possible  beam  reached  the  scattering  cell  beam  collector 
using  no  deflector  voltages. 
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CHAPTER  IV 
EXPERIMENTAL  RESULTS 
1.   PRELIMINARY  TESTS 
The  differential  scattering  cross  section  shown  in  equation 
(1)  of  Chapter  I  can  be  rewritten  in  terms  of  the  quantities 
actually  measured  as  follows: 

VG         2 

d<J<Q>  =   !.(p  P)t  an  cm  I1) 

lx   SC  ' 

where        I    =  current  at  MEM  detector  anode  in  amps 

G    =  MEM  detector  gain 

I.   =  incident  ion  beam  current  in  amps 

P    =  pressure  of  target  gas  in  Torr 
sc  ..  /- 

P    =  3.536  x  10   particles/cm  -Torr 
t    =  target  length  in  cm 
dQ  =  solid  angle 

To  measure  a  cross  section,  then,  a  target  gas  pressure  and  a 

target  length  are  established  in  the  system.   The  magnet  current 

and/or  the  detector  distance  are  varied  to  cover  the  range  of 

scattering  angles  to  be  observed,   These  two  quantities  also 

determine  the  value  of  the  multiplier  gain  and,  along  with  the 

angle  of  scatter,  the  solid  angle.   Then  at  each  setting  of  the 

magnet  current  and  the  detector  distance,  T   and  I.  are  observed 

D       1 

These  values  are  used  to  calculate  da  (9)  using  the  appropriate 
computer  program  listed  in  Appendix  IV. 

In  order  to  eliminate  any  spurious  signals  a  series  of  mea- 
surements was  undertaken  to  test  the  effect  of  the  variables  I., 

1 

P   and  t  upon  the  measured  value  of  do.      For  a  particular  angle 

of  scatter 

dcr  =  K  — — £ — -    ,    where  K    =    constant, 
i    sc 
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Now  for  constant  incident  ion  beam  current,  I.,  and  constant 
target  length,  t, 

*»■(¥)%•■  (2) 

The  experimental  data  plotted  in  Figure  21  verify  this  linear 
relation.   Notice  that  the  data  are  linear  but  does  not  extra- 
polate to  zero.   The  fact  that  I   does  not  go  to  zero  is  probably 

due  to  multiplier  background  noise.   However  this  error  is  not 

-3 
important  above  1.0  x  10   Torr  pressure  where  the  measurements 

were  usually  made.   Another  linear  relationship  which  must  exist 

■-  ■  (?  )■• 

for  constant  target  pressure  and  target  thickness.   Figure  22 
shows  a  plot  of  one  set  of  experimental  data  which  demonstrates 
this  characteristic.   Finally,  if  the  incident  ion  beam  current 
and  the  target  pressure  are  kept  constant  then  the  detector 
current  must  vary  linearly  with  the  target  thickness.   This  re- 
lation, which  is 

is  shown  experimentally  in  Figure  23.   These  three  curves  show 
that  our  detector  signal  has  the  proper  dependence  on  the  various 
parameters  and  behaves  like  a  scattered  particle  current.   For  the 

differential  scattering  cross  sections  measured  the  ion  current  at 

-7 
the  beam  collector  was  about  10   amps,  the  current  to  the  rear 

-9 

of  the  scattering  cell  was  about  10   amps,  the  ion  current  on  the 
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-8 
front  face  of  the  scattering  cell  was  about  10   amps,  and  the 

-8  -1? 

detector  signal  of  the  MEM  output  varied  from  10   amps  to  10 

amps.   With  no  target  gas  in  the  scattering  cell  and  essentially 
the  same  current  to  the  beam  collector  the  detector  current  drops 
to  about  10    amps.   Consequently  the  signal  due  to  scattered 
ions  was  at  least  100  times  greater  than  the  multiplier  dark 
current . 

2.   He+  +  He  SCATTERING  DATA 
Another  important  test  of  our  technique  is  the  comparison 
of  data  measured  by  our  method  to  data  obtained  by  more  con- 
ventional methods.   We  therefore  measured  the  elastic  differ- 
ential scattering  cross  sections  for  He  on  He  at  300  eV  from 
40   to  52   and  at  400  eV  from  42   to  54  ,  and  compared  our  data 
to  that  taken  by  Lorents  and  Aberth  ,  as  shown  in  Figures  24  and 
25.   The  gross  shape  of  our  data  and  Lorents'  data  is  essentially 
the  same,  but  we  do  not  observe  the  fine  structure  in  as  much 
detail  as  they  do.   This  is  due  to  two  effects,  both  tending  to 
increase  the  angular  resolution  (A0  —  1   for  our  data,  compared 
to  A0  —  0.5   for  the  Lorents  data),   One  effect  is  the  increase 
in  angular  spread  that  occurs  because  of  the  anomalously  large 
energy  spread  of  about  8,0  eV  in  the  incident  ion  beam  energy. 
The  other  effect  is  the  increase  in  angular  spread  brought  about 
by  the  spiraling  of  the  incident  ions  about  the  magnetic  axis. 
This  spiraling  will  allow  the  incident  ions  to  enter  the  target 
area  at  small  angles  with  respect  to  the  system  axis,  thereby 
increasing  the  angular  spread  the  detector  sees.   Nevertheless 
these  results  do  show  that  our  technique  can  give  reliable  cross 

section  measurements  on  an  absolute  basis. 
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3.   Li   +  He  SCATTERING  DATA 

Finally,  a  measurement  of  the  elastic  differential  scattering 

£i      + 
cross  section  of  Li   on  He  was  carried  out.   Since,  for  this 

system,  the  helium  target  atoms  have  a  smaller  mass  than  the 
incident  lithium  ions  there  are  two  possible  energies  at  which 
the  lithium  ion  may  scatter  at  any  allowable  laboratory  scat- 
tering angle.   There  also  is  a  maximum  angle  of  scatter  in  the 
laboratory  system  determined  by  the  relation 

°max  =  sin  (MHe/MLi>  =  ^-7°  <5> 

At  this  maximum  scattering  angle  the  laboratory  elastic  cross 

15 
section  becomes  infinite   .   This  particular  property  made  this 

system  suitable  for  a  check  of  the  accuracy  of  our  angular  de- 
termination and  system  alignment.   Since  our  experimental  appa- 
ratus does  not  detect  particles  scattered  at  a  particular  angle, 
but  rather  at  a  particular  vector  momentum  (recall  that  the  tra- 
jectory is  determined  by  the  p  x  B  forces),  this  peak  in  the 
laboratory  cross  section  can  be  seen  by  setting  the  detector  at 
a  convenient  distance  from  the  scattering  center  and  slowly  de- 
creasing the  magnet  current.   This  allows  us  to  measure  the  cross 
sections  continuously  from  the  high  energy  side  to  the  low  energy 

side  of  9    without  moving  the  detector  out  of  the  scattered 
max 

particle  beam.   Thus  we  are  able  to  determine  closely  the  magnet 
current  for  the  largest  cross  section  observed  and  compare  that 

magnet  current  value  with  current  value  for  0    calculated  from 

^  max 

the  theoretical  trajectories.   The  results  of  these  measurements 
are  shown  in  Figures  26  and  27  which  are  plots  of  the  laboratory 
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Figure   26.    da  (01    .     vs    Q,    .    for      Li      on  He,    E,    .    =   200   eV. 
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Figure   27.    da(°)lab  vs    Qlafe  for      Li      on  He,    E  =   300   eV 
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differential  scattering  cross  section  vs  the  angle  of  scatter  at 

200  eV  and  300  eV  respectively.   Notice  that  in  both  cases  the 

upper  branch,  corresponding  to  higher  scattering  energies,  suddenly 

increases  sharply  as  we  approach  0    (4l«7  )  and  then  decreases 

max 

along  the  lower  branch  as  the  angle  of  scatter  decreases.   The 
fact  that  the  cross  section  peak  appears  at  the  correct  angle 
illustrates  that  our  angular  determination  is  correct  and  our 
system  is  properly  aligned.   The  spread  in  the  data  is  about  a 
factor  of  2  and  is  believed  to  be  largely  due  to  error  in  the 
multiplication  factor  G  of  the  magnetic  electron  multiplier  and 
will  be  discussed  in  more  detail  below.   The  same  data  is  shown 
in  the  center  of  mass  system  in  Figures  28  and  29.   There  appears 
to  be  some  structure  in  both  curves  most  easily  seen  in  the  smooth 
data  in  Figure  30.   However,  from  this  data,  it  is  not  possible 
to  conclude  whether  this  structure  is  due  to  instrumental  effects 
or  represents  an  oscillation  in  the  cross  section  produced  by 

quantum  mechanical  interference  effects.   For  the  homonuclear 

+ 
case  (e.g.  He   +  He)  oscillations  produced  by  the  interference  of 

16 

gerade  and  ungerade  states  has  been  thoroughly  studied   ;   In 

+  17 

heteronuclear  systems  (e.g.  Li   +  He)  Lichten    has  proposed  that 

at  very  small  inter  nuclear  separations  a  similar  interference 

effect  may  occur  since  the  active  electron  might  find  itself  in 

a  field  resembling  that  of  a  homonuclear  molecule.   This,  of 

course,  would  cause  structure  in  the  differential  scattering  cross 

section.   On  the  other  hand,  there  are  reasons  to  believe  that 

instrument  effects  cause  the  oscillations.   In  the  80  eV  data,  the 

large  dip  in  the  cross  section  at  95   may  be  due  to  interaction 
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of  the  scattered  particle  beam  with  the  scattering  cell  exit  slit. 
This  problem  can  occur  for  particles  scattered  at  less  than  36 
in  the  lab.   The  dips  at  125°  in  the  80  eV  data  and  at  130°  in 
the  120  eV  data  correspond  to  particles  scattered  at  about  4l.7 
in  the  laboratory  system.   Consequently  it  is  possible  that  the 
dips  are  due  to  inaccuracy  in  the  calculated  solid  angle  very 
near  the  laboratory  cutoff  angle.   Notice  however  that  regardless 
of  the  structure  our  data  has  slightly  less  slope  and  is  above  the 

line  extrapolated  from  data  taken  at  smaller  angles  by  Lorents  and 

2 
Aberth  .   Our  data  indicates  that  the  center  of  mass  cross  section 

tends  to  level  out  as  the  center  of  mass  scattering  angle  increases 

Thus  there  is  essential  agreement  between  our  data  at  the  larger 

angles  and  that  of  Lorents  and  Aberth  at  the  lower  angles. 

k.       EXPERIMENTAL  ERRORS 

There  are  several  sources  of  error  in  the  measurements.   One 

source  of  error  which  cannot  be  entirely  eliminated  arises  from 

the  spiraling  of  the  incident  ions  about  the  focusing  magnet  axis. 

As  already  pointed  out  in  Section  2  above  this  introduces  an  error 

in  the  angular  measurement,  and  is  most  important  in  measurements 

where  the  cross  section  has  large  variations  as  a  function  of  the 

scattering  angle.   However,  the  effect  is  minimized  by  proper 

alignment  of  the  incident  beam  to  the  focusing  magnet  axis.   The 

gain  of  the  MEM,  which  is  the  scattered  ion  detector,  is  a  rather 

sensitive  function  of  the  magnetic  field  strength  of  the  focusing 

magnet.   Thus  any  small  error  in  setting  the  magnet  current  or  the 

detector  distance  would  result  in  actual  gains  slightly  different 

from  those  used  in  the  cross  section  calculations.   The  error  in 
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setting  the  detector  distance  and  the  magnet  current  would  also 
cause  differences  between  the  actually  observed  angle  of  scatter 
and  solid  angle  and  their  expected  values  from  which  the  cross 
sections  are  calculated.   The  detector  distance  can  be  set  to 

+  1  mm  and  the  magnet  current  can  be  regulated  to  +  0.05  amp. 

6   + 
For  the   Li   on  He  data  this  results  in  a  possible  error  in  the 

scattered  angle  of  less  than  0.05   and  in  the  solid  angle  of  less 

than  5%.   Another  source  of  error  is  in  the  measurement  of  the 

target  gas  pressure  in  the  scattering  cell.   Although  an  accurate 

determination  of  this  error  was  not  possible  it  is  apparent  from 

the  scattering  cell  geometry  that  the  actual  pressure  is  higher 

than  the  measured  pressure  used  in  the  cross  section  calculation. 

It  is  felt  that  the  primary  source  of  scatter  in  our  data  is  due 

to  the  multiplier  which  had  to  be  operated  at  potentials  different 

from  the  recommended  values  in  order  to  obtain  high  gain  in  the 

magnetic  field.   An  analysis  of  the  cross  sections  for   Li   on  He 

give  a  precision  which  is  on  the  order  of  +  20%  of  the  average 

value  of  a  set  of  measurements. 
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CHAPTER  V 
CALCULATION  OF  THE  INTERATOMIC 
POTENTIAL  OF  THE  Li+-  He  SYSTEM 
1.   INTRODUCTION 
The  estimation  of  the  interaction  energy  between  two  atoms 
as  a  function  of  their  separation  distance  is  of  primary  importance 
in  the  understanding  of  numerous  chemical  and  physical  problems. 
The  interaction  energy  can  be  found  in  several  ways;  however  the 
most  direct  method  is  to  calculate  it  from  elastic  scattering 
cross  sections.   In  this  chapter  the  potential  function  for  Li 

on  He  was  calculated  from  the  cross  sections  measured  by  our 

2 
technique  and  some  data  taken  by  Aberth  and  Lorents  .   The  pro- 
cedure used  to  calculate  the  potential  function  is  based  on 

3  18 

Firsov's   method  (outlined  by  Lane  and  Everhart   )  using  Gauss- 

Mehler  quaduatures  to  perform  the  necessary  integration.   The 

potential  function  was  then  fitted  to  a  Born-Mayer  form  and  com- 

19 
pared  to  the  potential  function  calculated  by  Zehr  and  Berry 

from  their  experimental  data, 

2.   THEORY 

Classical  scattering  theory,  valid  for  this  calculation   , 

relates  the  potential  function,  V(r),  to  the  center  of  mass  angle 

of  scatter,  0,  and  the  impact  parameter,  b,  by  the  expression 

0(E,b)=TT-2J   |[(1-  ^})r2-  b2l    dr  (1) 

r 
o 

where  E  is  the  center  of  mass  energy, 

r  is  the  separation  distance,  and 

r  is  the  distance  of  closest  approach, 
o  cir 
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The  necessary  and  sufficient  conditions  for  scattering  particles 

into  the  element  of  solid  angle  dQ   =  2rr  sinO  dO  is  that  the 

cm 

particles    be    incident    in    the   annular    ring   of   area    2nbdb,      There- 
fore  the    cross    section  cr(Q)    is    related   to   b  and   0   by 

2na(0)  sinO   dO   =    -    2nbdb 
or 

TT 

b2   =    2\   a (Q)    sinO   dO.  (2) 

0 

Since   the   functional    form   of  o"(0)    isn't    known   the    integration 
cannot    be    carried  out.      However    b   can   easily   be  found   from  a  plot 

of  cr(Q)    sin©   vs    0.      The   area   under    the   curve   from  rr    to   0   is    then 

2 
b  /2.      Thus    the   relation   between   0  and  b   can   be   obtained  directly 

from   the   scattering   data.      However    equation    (1)    cannot    be   used   to 

obtain  V(r)    directly,    so   the   following  method  will   be   used   to 

determine   the  potential   function    from  the   9   =   f(b)    result.      Let 


2        fi         V(t)n    2 


q  E 

OI  2       2 

V(r)    -   [1    -    qVr     ]  E  (3) 


and  by    knowing    q   and   r,    V(r)    is    determined.      Making    this    sub- 
stitution   in    (1)    gives 

0    =   n     -(    2b   [(^-)/dql  [q2-b2]      '    dq 


-  tt     -I    2b    (d  lr\   r/dq)(q2-b2)         dq 
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Now  using    the   fact    that 


=A      b(din    q2/dq)(q2-b2)      '    dq 
b 


leads    to 


-% 


i(b)     =    b)   [din(q2/r2)/dqil[q2-    b2]      '    dqx 


where   the    integration    variable   has    been    changed   to   q    .      Now 

2      2    2 
multiply    through  by   db/(b   -q    )      and   integrate   from   q  to    infinity. 


Thus 


oo  ,oo 


(b    -q^)2         J  J 

qv  ^    '  q  b 


[d^n(q2/r2)/dq]L]    ck^ 
2      2   3g      2      2   3g 

(b  -q  )   (q,  -b  ) 


(4) 


The    integral   on   the   right   can   be    transformed  using   Dirchlet's 

i     20    * 
formula         to  give 


0  ( b )  db 


1 


w  2 

'  d(qf) 


■C   - 


I  2      2.*g         2    \ 

(D  -q  )  J  2 

q    v  '  q^ 


—    d£„(-)dq 
U  r 


2      2.3g.    2    _2.3g 


q2    (b    -q    )':(q1-b    ) 


The   right   hand    integral    can  be   placed    in   a    standard   form   and 
integrated   to   give  rr  ,      Thus    equation    (4)    becomes 


Q  ( b )  db        _  n 
(b    -q    ) 


d(ql>  2      2  tt  2      2 

— —  d£n(q^/r2)dqi    =  |  in    l2/q2 


Thus 


r(q)    =    q   exp 


1    [       Q(b)db 


(5) 
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and  by  evaluating  this  equation  and  using  the  results  in  equation 

(3)  one  can  obtain  V(r).   The  evaluation  of  the  integral  in 

21 

equation    (5)    is    done   using   Gauss -Mehler    quadratures       ,    where 

n/2 


(    f(y)(l-y2)-%   dy   =  ^  £        f(cos   %!  tf) 


and  n  =  even  integer .   For  our  case  let  y  =  q/b,  then 


Thus 


Q(b)db  ^  =  f   Q(q/y)dy 
'q  (b2-q2)^  "  J0  y(l-yV 
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3.   PROCEDURE 

As  an  example  of  the  procedure  consider  the  calculation  of 

6   +  0 

V(r)  from  the  data  for   Li   on  He  at  80,0  eV .   The  cross  section 

is  plotted  in  Figure  30.   A  linear  fit  of  thisdata  was  used  in  the 
potential  calculation.   From  this  data  a  plot  of  a(0)sinO  is  con- 
structed.  This  plot  is  shown  in  Figure  31.   The  curve  has  to  be 

extended  to  180  .   The  area  under  the  curve  from  any  0  to  180 

2 
is  b  /2 .   The  result  of  this  integration  is  plotted  as  9(b)  vs  b 

and  is  shown  in  Figure  32.   From  this  graph  r(q)  is  evaluated. 
A  value  of  q  is  chosen  and  0(q/y)  is  found  for  each  term  of  the 
Gauss-Mehler  quadrature  formula.   Then  the  individual  f(y)  terms 
are  found  and  the  integral  is  evaluated.   From  this  one  obtains 
r(q)  and  finally  V(r)  using  equation  (3).   The  Gauss-Mehler  quad- 
ratures give  the  same  answer  to  four  significant  figures  for 

n  =  8,10,  or  12.   However  we  used  n  =  12  for  all  calculations, 

86 


100  120 

0rM  -    degrees 


'CM 


6    .+ 


Figure   31.      a(0)sin    0   vs   0  for      Li      on   He,    E        =   80.0   eV 
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Figure   32.      0(b)    vs   b  for      Li+   on   He,    Ecm   =   80.0   eV. 
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4.       RESULTS 

Using   the  procedure   outlined  above    the  potential   function 

was    calculated   from   our    data    at   E        =    80.0   eV   and   E        =    120.0    eV 

cm  cm 

and  from  Aberth  and  Lorents  data  at  E    =  109 . 1  eV  and  E   =  72.7  eV. 

cm  cm 

The  results  are  shown  in  Figure  33  where  the  potential  is  plotted 

against  the  internuclear  separation.   Notice  that  our  data  fits 

quite  well  with  that  of  Aberth  and  Lorents. 

The  plot  of  log   V(r)  vs  r  was  found  to  approximate  a 

straight  line,  thus  the  potential  function  could  be  represented 

by  the  Born -Mayer  form  as 

-4.4c        -8 
V(r)  =  259e      ,  r  in  10   cm. 

7      + 
Zehr    and   Berry   who  also  measured   the   scattering   of      Li      on   He 

found  their    data   to   fit   a    Born-Mayer    potential  which    is   expressed 
as  V(r)    =   370e"5olr 

This    is   plotted   on    the  same   graph  with   our   potential    (Figure   33). 

-8 
The   agreement    is    quite   good  down   to  about    r    =   0.3   x   10        cm, 

below  which   their    curve  rises   more   sharply.      This    discrepancy   at 

small    values    of   the    internuclear    distance    is    not    surprising    in 

view   of   the   fact    that   the   measurement   by   Zehr    and   Berry  was    carried 

out   by   a   difference  method  which   is   more   accurate   at   the   smaller 

angles    of   scattering   corresponding   to    larger    internuclear    distances 

The   agreement   with   the   data   of  Lorents    and  Aberth  appears    to  be 

excellent . 
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Figure  33.   Interatomic  Potential  for  Li   -He  as  a  Function  of 
Internuclear  Separation. 
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CHAPTER  VI 
CONCLUSIONS 

The  most  difficult  and  most  time  consuming  portion  of  the 
work  reported  here  was  designing  the  apparatus  and  demonstrating 
the  feasibility  of  using  the  focusing  properties  of  an  axially 
symmetric  magnetic  field  to  measure  differential  scattering  cross 
sections.   This  was  accomplished  by  showing  that  data  taken  by 
this  method  is  in  essential  agreement  with  that  measured  by  other 
means.   Therefore  this  method  with  its  much  larger  solid  angle 
will  allow  future  measurements  of  heretofore  too-small-to-be- 
observed  cross  sections. 

The  second  accomplishment  of  this  research  was  the  mea- 
surement of  the  differential  scattering  cross  section  of  Li 
on  He  at  large  angles,  thereby  extending  the  data  about  80   past 
previous  measurements.   These  measurements  also  demonstrated  the 
case  in  data  collection  that  this  technique  (with  its  built  in 
momentum  analysis)  offers  for  systems  with  heavy  incident  ions 
and  lighter  targets.   Finally  the  interatomic  potential  function 
at  distances  less  than  1A  for  the  Li   on  He  system  was  calculated 
from  our  data.   The  agreement  of  this  result  with  that  derived 
from  other  scattering  measurements  is  additional  evidence  for  the 
soundness  of  our  technique. 
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APPENDIX   I 
THE   EFFECT    OF  TARGET    SIZE   AND    ION    BEAM 
ENERGY   SPREAD  ON  THE  SOLID   ANGLE 
The   expression   for    the  angular    spread    (see   Section    5, 


Chapter    II)    was    /~ 


dzo  ■  \5S7         dE" 

dg    =  _ , ©R^ 

I        o 


sz  \  /3z 

EOZ'  "  EOR« 


}9    m>z< 


where      E      is    the    ion    beam   energy- 

dE    is    the  energy   spread   of   the  beam, 
dR'is    the   effective   target    radius, 
dz'is    the  effective   target    length. 

No   analytical   method  exists   for    calculating   the  four    partial 
derivatives.      They    are   computed  by   solving   a   number    of   trajectories 
where   E,Q,"R*,    and  Z1    are   varied  appropriately   and   the   corresponding 
change    in   the    intercept    distance,    Z    ,    is    noted.      The  value   of    dE 
offers    no  problem   since    it    is    the   energy   spread   of   the    incident 

ion   beam.      It    depends    only   on   the   source   being   used   and   can   be 

6     + 
measured   independently    (for    the     Li      source   dE  —   0„25  eV).      The 

quantities   dR  '    and   dz'    cannot   be  arbitrarily   chosen   because   they 
are   functions    of   the  trajectories   and   the   detector    geometry.      The 
maximum   size   of   the   interaction  region    is   determined  by    the   dia- 
meter   of   the    cell's    beam   entrance   aperture   and   the   cell's    exit 
slit    through  which   the   scattered    ions    depart   as    illustrated    in 
Figure   l8    (Section   5   of   Chapter    3).      However,    the  actual   region 
from  which   an    ion    can   scatter    and   enter    the   detector    is    smaller 
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than  this  and  depends  upon  the  trajectory  of  the  ion.   This  can 
be  seen  by  recalling  that  the  effective  detector  opening  is  de- 
termined by  AZ  ,  which  depends  on  0   (see  Appendix  II),  which  in 
turn  depends  on  the  particular  trajectory.   Thus  a  particle  can 

enter  the  detector  only  if  it  crosses  the  axis  within  Z   +  AZ  /2 

o  —   o 

and  at  the  proper  angle,  0  ,  regardless  of  where  it  originated 
in  the  target  gas.   Therefore  to  evaluate  dR T  or  dZ '  one  would 
have  to  evaluate  trajectories  for  many  different  values  of  dR1 
or  dZ*  until  the  trajectory  with  the  maximum  dR •  or  dZ ■  is  found. 
In  general  the  derivatives  RZ  /^R')  and  (bZ   /Bz*)  are  linear  for 

S>  Q  '  ^  O 

a  given  reaction.   Thus  the  change  in  dO  due  to  varying  values  of 
dR*  and  dZ'  is  linear  and  symmetric  about  its  value  for  dR ■  =  0 
and  dZ!  =0.   The  most  extreme  values  of  dR'  and  dZ *  an  ion  can 
have  and  still  be  detected  will  give  the  largest  positive  and 
negative  corrections  to  dO.   At  values  larger  than  these,  dO  is 
zero  and,  consequently,  dQ  is  zero.   Thus  it  is  apparent  that  the 
finite  target  size  corrections  are  such  that  the  average  dQ  is 
essentially  equal  to  the  value  of  dQ  obtained  for  scattering  from 
the  target  center.   It  is  suggested  that  a  detailed  study  of  the 
problem  be  made  on  an  analog  computer  by  tracing  the  trajectories 
backwards  from  the  detector  to  the  scattering  cell  and  defining 
dR'  and  dZ J  by  this  technique.   The  average  solid  angle  is  then 
obtained  by  integration  over  the  total  target  volume  contributing 
at  each  detector  and  magnet  setting. 

However,  one  may  obtain  an  estimate  of  the  maximum  effect  of 
finite  target  size  on  the  solid  angle  by  calculating  the  change 
in  AZ   due  to  changes  in  the  interaction  volume,  as  follows:  For 
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an  ion  to  be  detected,  it  must  recross  the  axis  within  a  length 

Z   +  AZ  /2  where  Z   is  the  distance  between  the  detector  and  the 
o  —   o  o 

scattering  center  and  AZ   is  determined  by  0   and  the  detector 

o  o 

geometry  (see  Appendix  II).   Now  let  the  detector  be  placed  at  Z 

so  that  it  will  detect  particles  scattered  from  the  center  of  the 

cell  at  an  angle  0.   The  angle  the  particles  make  with  the  axis 

as  they  recross  it  is  0  .   Ions  scattered  at  9  at  some  point  in 

the  target  other  than  the  center  will  recross  the  axis  at  Z'  with 

o 

an  angle  9r.   Since  AZ   decreases  as  9   increases  those  particles 
v  o  o  o  r 

with  9f  >  9  will  be  detected  only  if 
o    o 

AZ  AZ 

o  .  o 

Z r-  <  Z'    <  Z      +  -^- 

o     2  o  o     2 

and  there  will  be  no  increase  in  the  effective  solid  angle  from 

them.   However,  those  particles  with  9'  <  9  will  be  detected  when 

'       c  o    o 

AZf  AZ' 

z  -  -£-  <  z*  ^  z     +  -~ 

o     2      o     o     2 


and  the  effective  solid  angle  will  be  increased.   Therefore  we 

need  only  to  consider  those  trajectories  where  9 '  <  9  .   An 

'  J  o     o 

analysis  of  the  trajectories  of  ions  scattered  from  different 
places  in  the  target  show  that  those  scattered  on  the  ion  source 
side  of  the  center  and  those  scattered  off  the  axis  will  recross 
it  at  distances  less  than  those  scattered  from  the  center  and  that 
they  will  have  larger  9  values.  Consequently  they  cause  no  in- 
crease in  AZ  since  they  must  cross  the  axis  within  AZ  of  Z  to 
o  o     o 

be  detected,   Those  scattered  on  the  detector  side  of  the  center 
or  those  which  originate  off  the  axis  and  recross  it  before  leaving 
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the  scattering  cell  will  cross  at  distances  greater  than  those 

scattered  from  the  center  and  will  have  a  smaller  0   value.   Thus 

AZo  ° 
they  will  be  detected  if  they  have  a  Z1  ^  Z   +  — —  •   Therefore 
'  o     o     2 

the  actual  AZ   correction  is 
o 

AZ   +  AZ' 

A^C  0  0 

AZo  =  2 • 

Therefore,  for  a  particular  energy,  angle  of  scatter,  and  magnetic 

field  one  finds  (BO  /BZ  )  for  varying  dR ■  and  dz'  values  from 
v   o    o' 

trajectory  results,  then  gets 

BO 

e'  =  Q  +  (t^-)    az 
o    o    V3Z  '         o 

o 

which  gives  a  AZ } .   This  is  applied  to  give  the  proper  value  of 

c  c  6   + 

AZ  .   Table  II  shows  the  results  of  the  AZ   calculation  for   Li 
o  o 

on  N  at  100  eV,   Notice  that  varying  dR '  gives  the  larger  value 

for  (BO  /Bz  )  and  therefore  the  smallest  value  of  0'.   This  was 
v   o    o'  o 

used  to  calculate  AZ*  since  this  was  the  largest  correction.   The 

o 

largest    error    is    about    8%  and   as    the   angle    increases    the   error 
decreases  '. 

The   energy   correction    is    straightforward.      The   derivative 
(Bz   /BE)    is    evaluated   from   trajectory    data.      The   energy    spread 
dE   is   determined   separately   for    the  particular    ion   source   and 
the   correction    (BZ  /BE)dE    is    calculated.      This    correction    is 
about    6%. 
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APPENDIX    II 

CALCULATION   OF  THE   INTERCEPT    SPREAD 

The   fact    that    the   detector    apertures    are   of   a    finite   size 

allows    all    ions    whose    intercept    values    are   Z      +    AZ  /2   to   enter 

r  o  —   o 

the  detector  and  be  collected.   This  allowable  intercept  spread, 
AZ  ,  depends  on  the  detector  aperture  geometry  and  on  the  angle 

Q   between  the  ion  path  as  it  enters  the  detector  apertures  and 

o  r  r 

the  center  line.   This  is  readily  seen  in  Figure  A2-1.   In 

general,  then,  as  0   increases,  AZ   decreases. 

o  '         o 

The   specific   relationships    between    AZ      and   0      can   be   found 

c  o       o 

as  follows : 

Consider  an  incident  beam  of  ions  whose  outer  edge  is 
determined  by  the  outer  aperture  and  whose  inner  edge  is  deter- 
mined by  the  multiplier  grid  opening.   AZ   for  this  case  can  be 
calculated  from  the  geometry  shown  below. 


-b— 

— b- 

, 

•Vs. 

^S 

f\^ 

^6o 
1 

_► 

AZ( 

3 

Now 
but 

and 


AZ   =  (-4+a)  cotan  0 

0  v      '         ( 

1  =    s  cotan  0 

o 

s+f  _  .    n 
-— -  =  tan  0 
2b         o 
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MULTIPLIER 
CATHODE 


MULTIPLIER 

ENTRANCE 

GRID 


ELECTRON 
MULTIPLIER 


b  =  1.27cm 
f  =  0.95cm 
h  =  0.32cm 


TA-370S22-I 


Figure  A2-1.   Detector  Aperture  Geometry. 
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or  s  =  2b  tan  0  ~f 

o 

giving         I    -    2b-f  cotan  0 
y     *  o 

which  gives   AZ   =  (f+a)  cotan  0  -2b.  (1) 

o  o  , 

The  maximum  0   a  particle  can  have  is  that  angle  defined  by  a  line 

from  the  outer  aperture  edge  to  the  farthest  edge  of  the  multiplier 

grid.   For  this  trajectory  AZ   =0.   Thus 

AZ   =  0  -  (f+a)  cotan  0  -2b 
o  '  o 

cotan  0   =  2b/ (f+a) 
o 

i   2b 
0   =  cotan-1  (!^-) 
o  vf+a; 

putting  the  values  shown  in  Figure  A2-1  gives  0  .      .  =  33-2  . 
,  .  o(max  I ) 

As  we  decrease  0  -equation  (1)  holds  until  the  angle  of  the 
incoming  ion  is  the  same  as  the  angle  formed  by  the  edge  of  the 
outer  aperture  and  edge  of  the  shield  aperture.   This  angle  is 
found  from  the  relation 

tan  0   =   , 
o    b 

_        -1  ,f-h. 

0   =  tan    f-r— ) 
o  v  b  • 

and  putting    the   values 

0         .  =0  =   26,6°, 

o(min    I)  o(max   II) 

This  is  the  smallest  0   for  which  equation  (1)  determines  AZ  . 

o  v  '  o 

Now  consider  an  incident  beam  of  ions  with  0   =26,6   whose 

o 

outer  edge  is  determined  by  the  shield  aperture  and  whose  inner 
edge  is  determined  by  the  multiplier  grid  opening.   This  is  shown 
in  the  following  figure. 
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Again 


and 


giving 


AZ   =  -  i  +  '\a  cotan  0 


o 

I  = 

s    cotan    0 
o 

h+s 
b 

tan    0 
o 

s    = 

b   tan   0      -h 
o 

AZ   = 


(h+a)  cotan  0   -  b, 
v    '  o 


(2) 


As  9   decreases  the  limit  on  the  inner  edge  of  the  incident  ion 

o  9 

beam  will    finally    be   determined   by    the   line   from   the   edge   of   the 
shield  aperture   to   the   nearest    edge   of   the   multiplier    grid.      At 

this   time 

n  a-h 

tan    0      =   — -— 
o  b 


or  0     .     .  =    Q  =    17.3    . 

o(mm    II)  o(max   III)  J 


This  is  the  lower  limit  on  the  use  of  equation  (2). 

The  final  case  occurs  when  the  incident  beam  of  ions  is 
limited  only  by  the  shield  aperture.   The  following  sketch 
illustrates  these  limits. 
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Here  AZ      =    2h    cotan    0    .  (3) 

o  o  v     ' 

Thus,    the   three  possible  ways    a   beam   of    incident    ions    can  be 
limited  by   our    detector    geometry  along   with  the    applicable   ex- 
pression   for    AZ      are: 
r  o 

I.    For      33.2°  >    0     >    26.6° 

o 


AZ      =    (f+a)    cotan    0    -2b; 
o  o 


II.    For       26.6°  >    0     ^    17.3° 

o 

AZ      =    (h+a)    cotan    0    -b; 
o         v  '  o      ' 


III.    For      17.3°  >    0     >    0° 

o 


AZ      =    2h   cotan   0 
o  o 


Figure  A2-2    shows    a   plot    of   AZ      vs    0      over    the   appropriate 


angular    range. 
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Figure  A2-2.   Intercept  Spread  vs  Intercept  Angle 
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APPENDIX  III 
CALCULATION  OF  TARGET  THICKNESS 
The  thickness,  t,  of  the  target  depends  on  the  scattering 
angle,  0,  the  slit  openings,  which  can  be  controlled  from  outside 
the  vacuum  system,  and  the  angles  which  the  sides  of  the  scattering 
cell  make  at  the  exit  slit. 


a  =  2.5k  cm 
b  =  2.00  cm 
h  =  1.50  cm 


From  the  above  sketch  of  the  scattered  particle  beam  exit, 
the  four  points  (numbered  1,2,3,4,)  that  restrict  the  scattered 
beam  can  be  seen.   For  particles  scattered  at  angles  less  than 
36  ,  points  2  and  4  limit  the  scattering  length.   For  particles 
scattered  at  angles  equal  to  or  greater  than  49  ,  points  1  and  3 
limit  the  scattering  length. 

Now  consider  t  for  particles  scattered  at  angles,  0,  less 
than  36  .   Then  t  must  be  the  distance  along  the  cell  axis  between 
two  parallel  lines  each  making  an  angle  0  with  the  axis,  where  one 
of  the  parallel  lines  passes  through  point  2;  the  other  through 
point  4,  as  shown  on  the  following  sketch. 
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To    find   t    analytically,    construct    a    third    line  parallel    to  the 

first    two   and   require   that    it   pass    through  point    1. 

Then  t    =    s    -  i    and  I    can  be  found   from  the    law   of   sines 


I 


sin(36    -0) 


sin    0 


or 


I    =   a    sin (36    -0)/sin   0 


and 


t    =    s    -    a    sin(36    -0)/sin    0 


(1) 


Now   consider    t    for    particles  scattered   at    angles    0   equal   to 

y-O  I         o 

or    greater    than    3°      but    less    than  or    equal   to   £+9    •      This    case    is 
shown    in    the    sketch   below. 
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Here  t  is  the  distance  along  the  scattering  cell  axis  between  two 
parallel  lines  each  making  an  angle  0  with  the  axis;  one  of  the 
lines  passing  through  point  1,  the  other  through  point  4.   Then 
it  is  obvious  that  for  this  case 

t  =  S.  (2) 

Finally,  consider  t  for  particles  scattered  at  angles  greater 
than  49  •   For  this  case  t  must  be  the  distance  along  the  cell 
axis  between  two  parallel  lines  each  making  an  angle  0  with  the 
axis;  one  of  the  lines  passing  through  point  1,  the  other  line 
passing  through  point  3,  as  shown  below. 


It  can  be  seen  from  the  sketch  that  t  =  S  -  &. 
Again  using  the  law  of  sines  one  has 


i 


Sin(0-49O)    sinO 


giving 


I    =   b  Sin(0-49  )/SinO 

t  =  S  -b  Sin(0-49°)/SinO. 


In  summary:     0  <  36°:  t  =  S  -a  Sin (3o°-0)/SinO 
36°  <  0  <  49°:  t  =  S 


49°  *  0 


t  =  S  -b  Sin(0-49  )/SinO 


(3) 
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A  plot  of  i    vs  Q  is  shown  in  Figure  A3-1-   From  this  relation 
the  scattering  length,  t,  can  be  obtained.   Figure  A3-2  is  a  plot 
of  t  vs  0  for  a  scattering  cell  opening  of  0.3  centimeters.   Note 
that  it  depends  on  two  things:  (1)  the  length  of  the  scattering 
cell  exit  slit  and,  (2)  the  angle  of  scatter  being  observed. 
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APPENDIX  IV 
LISTING  OF  COMPUTER  PROGRAMS 
The  computer  programs  used  to  calculate  the  solid  angles  and 
the  cross  sections  are  listed  in  this  appendix,   All  the  programs 
were  written  in  Fortran  VI  and  were  run  on  the  IBM  36O-67  computer 
The  programs  are: 

1.  SOLANG:   This  is  used  to  integrate  the  trajectory 
equation  and  to  calculate  the  solid  angle.   This 
program  was  written  by  N.R.A.  Smyth. 

2,  CSVSZ:   This  program  calculates  the  differential 
scattering  cross  section  from  the  raw  data  where 
the  observed  angle  is  varied  by  moving  the  detector 
while  keeping  the  magnet  current  constant. 

3-   CSV SI :   This  program  calculates  the  differential 
scattering  cross  section  from  the  raw  data  where 
the  observed  angle  of  scatter  is  varied  by  changing 
the  magnet  current  while  keeping  the  detector  fixed. 
4.   PSIGMA:   This  program  calculates  the  differential 
scattering  cross  section  from  the  raw  data  for  the 
case  where  the  target  has  less  mass  than  the  pro- 
jectile. 
All  three  cross  section  programs  give  both  the  laboratory  system 
and  center  of  mass  system  values  for  the  cross  section  and 
scattering  angle. 
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APPENDIX  V 

ION  BEAM  OPTICS 

22 

Many  good  references    on  the  design  and  use  of  electro- 
static and  magnetic  lenses  to  focus  ion  beam  already  exist  and 
are  well  known  by  workers  in  the  field.   Therefore,  only  the 
space  charge  effect  and  the  focusing  properties  of  an  einzel  lens 
will  be  discussed. 
Space  Charge 

The  coulomb  interaction  among  the  charged  particles  in  an 
ion  beam  causes  the  beam  to  spread  as  it  moves  away  from  the 
source.   This  spreading  due  to  the  space  charge  sets  an  upper 
limit  on  the  amount  of  beam  that  can  be  transported  between  two 
apertures.   The  addition  of  converging  lenses  along  the  path  of 
the  ion  beam  will  help  reduce  the  space  charge  spreading  and  help 
deliver  more  beam  to  the  target  area. 

The  maximum  current  that  can  be  transported  a  distance  x 
starting  with  a  radius  y   and  ending  with  a  radius  y   can  be 

C  3l 

calculated  from  the  following  assumptions: 

(1)  The  radial  current  density  is  constant  across  any 
beam  cross  section. 

(2)  The  current  density  is  zero  outside  the  outer  ray. 

(3)  The  total  current  is  constant. 

(4)  The  axial  velocity  of  the  beam  is  constant. 
The  field  E  at  the  surface  of  the  beam  can  be  calculated  from 
Gauss'  Theorem.   Thus 

jpdv  =  fe  E   da 
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where  P       =    space    charge   density 

da    =    surface    element 
dv    =    volume   element. 

Now    integrating   at    y      gives 

2 

P  H  y       Ax    =    €     E       2TTy  Ax 

c  y 

since   all    the   flux   goes    out    through   a    cylindrical    surface    at 

radius    y.      Thus 

2  . 
E     =  Py  /2ye 
y  c 

2 
The    current    1    through   a    cross    section  TTy      can   be   found   from   the 

c 

linear    velocity    of   the    ions,    u    c    Therefore 

A,      2 

I  Ay     =  pu 
c  x 


or 


2 
P    =   i/ttu  y 

X    c 


2 

Now                            _                d   y  /„    2  . 

qE      =   m  y-   =  q(py    /2ye  ) 

y              dt  C 


and    substituting   for    p    and   solving   gives 


^   =   -SL-    (-L-)    =   b2/2y 
2        2mu      vney;  '  ^y 

dt  x 


where 


Now 


2 
b      =   ql/mu  ne    -    constant. 
^  x 


J 


fiK*  «•  ■  £  f  *)-■ 


integrating       /,    \2 
/dy  ] 

dt  / 


dZ^     =    b2    log   y+c,  (1) 
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The  initial  radius  was  y   and  at  that  time  all  ions  were  traveling 
parallel  to  the  x  axis,  thus 

dyc  2 

— —  =  0  giving  c  =  -b   log  y  . 

dt      9-9  s  c 

Substituting  into  equation  (1)  gives 


x 
dt         ~v-~»   "  »  c 


-2-  -   b(log   y/y    )2  which    is    the   differential   equation 


for    ions    in   the   outer    ray.      To    solve   this    equation    let 

h 
s  =  (log  y/yc) 

or  2 

y    =    yceS       . 


Then 


2 
dy  ,dywds,         ,  _         s      ds 

dt         <dt><dt>    =   bs    =   yc   2s    e        — 


giving  2 

s  ds 

b   =   2y      e  — 

c  dt 


Now  t  s 

i  b       --         C         s      ds 

dt    =   \      e 


or 


2y  )  dt 

t         C  Jo 

c 


<*-*<:>  = )   e     dT  =  D(s)-D(sc) 


2y 

o 

23 
where   D(s)    is    the    Dawson   function,    which   has    been   tabulated      .       But 

t-t      =    (x-x    )/u    ,    which 
c         v         c'      x 


gives 


x-x 

(- —)    b   =    D(s)  . 

v2y    u    '  v     ' 

ex 
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2u 

Thus  (x-x  )/y   =  — ~   D(s). 

v    c'   c     b     v  ' 

x, 
2qV  « 

Now  using  u   =  (— -1— )   and  also  substituting  for  b  gives 

9     x    v  m  '  "       " 


4(™)   (£) 


(x_Xc)/yc  =  o/p  I/""    D(s)« 

Therefore  the  maximum  ion  beam  current  that  can  travel  a  distance 
starting  at  a  radius  y  and  ending  with  a  radius  y,  whose  energy 
is  qV,  for  ions  having  a  charge  q  and  mass  m  is 

y     2        2 

I  =    I6ne  (q/2m)^  V3/2(^r)   (D(s))  .  (2) 

c 

For  our  experimental  apparatus 

x-x   =  75.0  cm 
c 

y   =  0,25  cm 
c 

y  =  0.25  cm 

2 
s   =  log  y  /y  =  1  and  D(s)  =  1.2463,  giving 

-9   V3/2 
I  =  5.38  x  10  v  V- 

m^ 
where  I  is  in  amps,  V  is  in  volts  and  m  is  in  amu. 

A  plot  of  I  vs  V  is  shown  in  Figure  A5-1  for  several  different 
ions . 
Binzel  Lens 

In  order  to  reduce  the  space  charge  spreading  electrostatic 
lenses  were  placed  at  the  source  and  also  about  half  way  between 
the  mass  analyser  and  the  scattering  cell.   The  three  diaphram 
einzel  lens  shown  in  Figure  A5-2  was  chosen  for  both  lenses. 
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Figure  A5-2.   Einzel  Lens  Equipotential  Lines. 

In  this  lens  the  two  outer  diaphrams  are  set  at  ground  potential 
and  the  potential  on  the  inner  diaphram  is  varied  to  give  the 
maximum  current  at  the  scattering  cell.   Since  the  voltages  on 
both  sides  of  the  lens  are  the  same  the  refractive  indices  are 
equal  on  both  sides  and  the  three  diaphrams  act  as  a  single  lens 

with  one  focal  length.   In  fact  the  focal  length  was  found  to  bt 

24 
independent  of  the  applied  potential.   Chancon    investigated  this 

lens  by  ray  tracing  and  found  the  following  empirical  formula  for 

the  focal  length 

i.  6.5  /f  tan'1(f» 

s       I  ,    s  .   -l.s 


1+Rtan   <R> 


where   s  =  lens  separation 

R  =  center  aperture  radius 
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The   distance  between  principal   point   P    and  mid-plane  M    is    given   by 

Rs 


MP    =    0.06 


R    +T 
o 


where  R   =  outer  aperture  radius 

T  =  thickness  Of  the  diaphrams. 
The  focal  length  relation  holds  as  long  as  T  <  — jr  S.   A  plot  of 
the  focal  length  relation  is  shown  in  Figure  A5-2.   For  the 
lenses  used  in  our  apparatus 


S  =  1.4  cm 

R  =   0.7  cm 

which  gives   f  =  6.3  cm 
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APPENDIX  VI 
ELECTRON  IMPACT  ION  SOURCE 

At  the  beginning  of  this  experiment  it  was  realized  that 
several  ion  sources  compatible  with  the  experimental  apparatus 
would  be  required.   Therefore  some  extra  effort  was  put  into 
building  an  electron  impact  ion  source  and  determining  its  operating 
characteristics.   Although  the  majority  of  the  data  presented  in 
this  thesis  were  obtained  using  a  different  ion  source  it  is 
planned  to  use  plasma  type  ion  sources  in  future  measurements  and 
a  summary  of  our  developmental  work  is  therefore  included  in  this 
appendix . 

In  this  ion  source,  shown  schematically  in  Figure  A6-1,  the 
ions  are  produced  by  electron  bombardment  of  the  gas  filling  the 
source  can.   A  plasma  is  formed  in  the^-region  between  the  grid 
and  the  filament  and  the  ions  are  extracted  axially  from  it  by 
the  application  of  external  electrostatic  fields. 
Physical  Description 

The  ion  source  consisted  of  a  5.715  centimeter  diameter  in- 
verted glass  press  (see  Figure  A6-1)  made  of  number  3320  glass, 
with  four,  in  line,  equally  spaced,  1/16  inch  diameter  pure  tung- 
sten rods.   The  two  outer  rods  support  a  15  mil  tungsten  wire 
helical  grid.   The  two  inner  rods  support  a  15  mil  tungsten  wire 
filament.   Both  grid  and  filament  are  4*8l  centimeters  long  and 
4.8l  centimeters  above  the  flat  side  of  the  glass  press.   In  the 
original  configuration  the  glass  press  became  coated  with  a  layer 
of  tungsten  after  a  very  short  period  of  operation.   In  order  to 


140 


Deflector  Shims 


Insulators 


Supports 


Lavite  Ring 


Einzel  lens 


Extraction  Aperture 


Tungsten  Rods 


^^         ffT '^-^ — Glass  Press 
^■-^        J  ^—Source  F 

P^77777l 


lange 


Kovars 


Gas  Inlet  m         YA 

Ion  Gauge  Fitting 


0  I 

L   i   i   i    I 


2  3 

i  I 


SCALE    INCHES 


Figure  A6-1.      Schematic   of  Electron    Impact    Source. 
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prevent  the  filament  and  grid  from  shorting,  the  glass  press  was 
inverted,  letting  the  tungsten  support  rods  sit  in  the  well  as 
shown  in  the  figure.   With  this  arrangement  no  shorting  of  the 
source  elements  was  experienced. 

Source  gases  were  passed  through  a  liquid  nitrogen  trap  to 
remove  water  vapor  before  they  entered  the  source.   A  precision 
gas  metering  valve  placed  in  the  inlet  pipe  was  used  to  regulate 
the  gas  flow  into  the  source.   The  source  pressure  was  measured 
with  a  Bayar d-Alper t  type  ion  gauge  attached  to  the  source  by 
means  of  a  3/4  inch  inside  diameter  "quick-connect"  coupling 
attached  to  the  source  flange.   This  arrangement  gave  positive 
control  and  accurate  measurement  of  the  ion  source  pressure. 

The  source  can  was  made  of  number  304  stainless  steel  and 
insulated  from  its  base  plate  with  a  1.27  centimeter  by  1.27 
centimeter  lavite  ring.   The  beam  extraction  hole  was  varied  in 
size  from  0.5  centimeter  to  0.1  centimeter.   An  einzel  lens  con- 
sisting of  three  parallel  stainless  steel  plates  was  attached 
directly  in  front  of  the  ion  source.   The  holes  in  the  lenses  were 
1.4  centimeters  in  diameter  and  the  lens  elements  were  1.27  cen- 
timeters apart.   These  lens  dimensions,  calculated  using  the  lens 
formulas  given  in  Appendix  V,  were  found  to  give  the  highest  beam 
currents  inside  the  scattering  cell.   Two  pairs  of  steering  elec- 
trodes are  mounted  external  to  the  outermost  focusing  lens  to 
assist  in  beam  alignment. 

The  ion  source  is  wired  as  shown  in  Figure  A6-2,   The  fila- 
ment power  is  supplied  through  an  isolation  transformer  which 
allowed  the  filament  and  the  grid  to  float  on  top  of  the  extraction 

142 


o> 


^M 


C     4 


6 

H 

o 

•H 

Q 

o> 

c 

•rl 
M 
•H 

01 

o 

VI 

0 
§ 

M 

o 

<fl 

& 


W 

«-» 
U 

01 

rH 


I 

< 
0) 


•H 

u. 


143 


electrode  (can)  potential.   The  highest  potential  is  that  of  the 
grid  and  is  the  approximate  beam  energy.   Energy  analysis  has 
shown  the  energy  spread  to  be  +  4-0  eV  about  the  gr id  potential, 
however  other  workers  have  found  the  energy  spread  to  be  +  0.5  eV < 
A  filament  current  of  approximately  11  amps  was  found  to  produce 
adequate  electron  emission  of  about  100  ma.   This  corresponds  to 
a  measured  filament  temperature  of  about  38OO  K  and  a  filament 
lifetime  of  over  100  hours. 
Operating  Characteristics 

The  grid-current/grid-voltage  relationship,  shown  in  Figure 
A6-3  for  H   gas,  are  characteristic  of  these  of  all  gases  used. 
The  source  operating  area  corresponds  to  the  grid  voltages  above 
the  grid  current  plateau.   Although  there  is  some  measurable  ion 
beam  with  a  zero  grid  voltage,  there  appears  to  be  a  minimum  grid 
voltage  below  which  there  is  no  actual  discharge.   These  minima 
were  found  to  be,  for  the  different  gases, 

GAS  MIN.  GRID  VOLTAGE 
H  70  Volts 

He  90  Volts 

N2  36  Volts 

-3 
at  a  source  pressure  gauge  reading  of  3.0  x  10    Torr .   Any  in- 
crease above  these  voltages  caused  a  small,  1%   at  most,  decrease 
in  total  beam  intensity. 

The  dependence  of  ion  beam  intensity  on  source  gas  pressure 
is  shown  in  Figures  A6-4  and  A6-5  for  hydrogen  and  nitrogen  re- 
spectively.  As  can  be  seen  from  the  figures,  the  maximum  ion 
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O    -   data  taken  by  decreasing  voltoge 
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Figure  A6-3.      Electron   Impact   Source   Characteristics. 
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beam  intensity  occurs  at  source  ion  gauge  readings  between 

-3  -3 

2  x  10    Torr  and  3  x  10   Torr,   Below  this  pressure  range  the 

ion  beam  intensity  is  quite  small  since  there  was  no  longer  a 

discharge  in  the  source.   With  the  ion  source  operating  properly, 

-8 
beam  intensities  on  the  order  of  3.0  x  10    amps  were  obtained  in 

the  scattering  cell. 

Theory  of  Operation 

The  discharge  in  the  source  is  a  Townsend  discharge,  which 
is  non-self-sustaining.   There  are  two  primary  sources  of  ionizing 
electrons;  one  is  the  tungsten  cathode  which  produces  electrons 
through  thermionic  emission,  the  other  is  neutral  atoms  or  mole- 
cules from  which  electrons  are  liberated  in  ionizing  collisions. 
However,  before  these  "secondary"  electrons  can  ionize  in  turn 
they  must  be  accelerated  by  the  electric  field  until  they  obtain 
a  kinetic  energy  equal  to  or  greater  than  the  ionization  potential 
of  the  neutral  gas.   The  effect  of  these  secondary  electrons  can 
be  seen  in  Figure  A6-3.   There  the  large  increase  in  grid  current 
and  beam  current  as  the  grid  voltage  is  increased  from  60  volts 
to  90  volts  is  explained  by  the  fact  that  above  about  60  volts 
the  secondary  electrons  have  enough  energy  to  ionize  the  neutral 
gas  through  collisions. 

The  beam  current  is  also  pressure  dependent  as  is  seen  in 
Figures  A6-4  and  A6-5.   This  peak  in  the  beam  current  vs  pressure 
curve  can  be  explained  in  the  following  way.   An  electron  with  a 
small  energy  in  an  electric  field  E  begins  a  mean  free  path,  L. 
To  be  capable  of  ionizing  at  its  next  collision  it  must  gain  energy 
from  the  electric  field  after  traveling  a  distance  d  such  that 
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v 

qdE  ^    q      or    d  ^    — ■ 

where   V    is    the    ionization   potential    of    the   neutral   gas.      The 
probability    that    the   electron  will    travel   a    distance   d   without 
a    collision    is 

P(d)    =   e"d/L    . 

This  is  the  number  of  ions  produced  by  an  electron  in  a  mean  free 
path.   Since  1  cm  of  actual  path  contains  1/L  mean  free  "ionizing" 
paths j  the  probability  of  a  free  path  longer  than  L  per  cm  of  path 
in  the  field  direction  is  equal  to  the  number  of  ionizing  col- 
lisions per  cm  of  path  in  the  field  direction.   Therefore  the 
number  of  ions  produced  by  an  electron  per  centimeter  is 

a    =  p(d)/L  =  i  e~d/L 

1   "el 

a    =  —   e 

The  mean   free  path   depends    on   the  pressure   and    can   be   written   as 

—  =  Ap ,    where  A   is    a    constant    and  p    is    the  pressure.      Then 

L  JL 

Ep 
Ae 
a   =  p  ,    where   B  =  AV. 

From   this    expression,    the  pressure   at   which  a    is    a   maximum, 
making   the   beam   current   a   maximum,    can   be   found. 

.        -B/Ep     ,  .    B       -B/Ep     . 

da  =   A   e  r   dp    -   p  A  -  e  *dp 

a 


dot        ..         B      , 

d?=    <X    -    E   P) 
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At  maximum    da  .  _J     >    „ 

-7—  =  0  making  (1  -  B/E  p  )  =  0, 

dp  9  rm' 

Thus  p   =  E/B 

is  the  pressure  at  which  maximum  beam  current  (ion  formation) 
occurs.   Notice  that  this  depends  on  the  gas,  its  ionizing 
potential,  and  the  electric  field  the  electron  is  in.   This 
demonstrates  that  there  should  be  a  maximum  in  the  pressure  vs 
beam  current  data.   However,  our  approach  has  been  too  sim- 
plified to  allow  us  to  calculate  p   from  tabulated  values  of  B, 
r  ^m 

For  example  ionization  by  photons  emitted  from  the  tungsten 

filament  has  been  ignored.   However,  from  the  simple  model  it 

is  seen  that  in  order  to  obtain  the  most  intense  beams,  the 

source  should  be  operated  at  or  a  little  above  p   and  with  a 

*m 

grid  voltage  at  or  above  that  required  to  produce  ionizing 
secondary  ions. 
Proton  Production 

Since  a  high  intensity  proton  beam  is  desirable  for  the 
study  of  some  of  the  proposed  reactions  an  attempt  was  made  to 
enhance  the  proton  output  of  the  electron  impact  ion  source  by 

introducing  water  vapor  into  the  discharge  as  was  suggested  by 

25 
Dawson  and  Tickner   .   The  set  of  reactions 


H3  +  H2°  *-*  H3°+  +  H2 


H  0+  +  e"— ■*  H20  +  H 
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if  actually  occur ing,  would  enhance  the  amount  of  atomic  hydrogen 

+ 
and  hence  the  proton  output  while  reducing  the  H„  output.   It  was 

found  that  the  proton  output  was  increased  with  a  corresponding 

+       + 
reduction  in  the  H   and  H„  outputs  as  water  vapor  was  added  to 

+ 
the  discharge „   However  the  total  H  beam  intensity  itself  de- 
creased by  a  factor  of  100.   It  was  also  found  that  significant 
changes  in  the  ratio  of  H  to  H„  and  H_  yields  occur ed  only  after 
the  discharge  was  20%  to  k0%   water  vapor.   It  was  felt  that  this 
was  an  unacceptable  load  on  the  vacuum  pumps  and  this  approach 
was  therefore  not  pursued  further. 
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